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Abstract 

The current status of the physics of massive neutrinos is reviewed with a 



00 . 

. forward-looking emphasis. The article begins with the general phenomenology 

m ■ 

■ of neutrino oscillations in vacuum and matter and documents the experimental 

' evidence for oscillations of solar, reactor, atmospheric and accelerator neutri- 

^ \ nos. Both active and sterile oscillation possibilities are considered. The impact 

1^ . of cosmology (BBN, CMB, leptogenesis) and astrophysics (supernovae, highest 



energy cosmic rays) on neutrino observables and vice versa, is evaluated. The 
predictions of grand unified, radiative and other models of neutrino mass are 
discussed. Ways of determining the unknown parameters of three- neutrino os- 
cillations are assessed, taking into account eight-fold degeneracies in parameters 
that yield the same oscillation probabilities, as well as ways to determine the 
absolute neutrino mass scale (from beta-decay, neutrinoless double-beta decay, 
large scale structure and Z-bursts). Critical unknowns at present are the am- 
plitude of Ue oscillations and the hierarchy of the neutrino mass spectrum; 
the detection of CP violation in the neutrino sector depends on these and on an 
unknown phase. The estimated neutrino parameter sensitivities at future facil- 
ities (reactors, superbeams, neutrino factories) are given. The overall agenda 
of a future neutrino physics program to construct a bottom-up understanding 
of the lepton sector is presented. 
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1 Introduction 



After decades of immense experimental and theoretical effort, major breakthroughs 
in our understanding of the properties of neutrinos have occurred recently. Since 
the time that the neutrino was first proposed by Pauli PP and placed on a concrete 
theoretical foundation by Fermi [2] , the question of whether neutrinos are massless or 
massive has persisted^. The standard technique of probing neutrino masses through 
studies of the endpoints of decay spectra succeeded only in placing increasingly more 
restrictive upper limits on neutrino masses, presently down to about 2 eV in the case 
of neutrinos emitted in the beta decays of tritium jnHEj- 

A more sensitive measure of small neutrino masses was known since 1968 from the 
proposal of Gribov and Pontecorvo that a neutrino of a given initial flavor could inter- 
change its identity with other flavors |7j, with a probability that is dependent on the 
distance from the location of the source |H]. Early searches for evidence of neutrino 
oscillations found only upper bounds on the oscillation probabilities. A long-standing 
deficit of solar neutrinos in the '^''Cl radiochemical experiment of Davis [HI, compared 
to Standard Solar Model (SSM) expectations [TUillll, was attributed to oscillations, 
but this interpretation remained unproven until recently. The observation of an elec- 
tron to muon ratio in the atmospheric neutrino events in the Kamiokande [12| and 
1MB ^] experiments of about a factor of 2 above expectations was interpreted as 
evidence for oscillations with neutrino mass-squared difference 5m? ~ 10^^ eV^ and 
near maximal neutrino mixing j^. However, due to the prevailing theoretical preju- 
dice that neutrino mixings would be small, based on the fact that quark mixings are 
small, this interpretation of the atmospheric neutrino data did not receive widespread 
acceptance. 

The definitive evidence of atmospheric neutrino oscillations came ten years later 
from the Super-Kamiokande (SuperK) experiment ^3]. With the ability to measure 
the zenith angular and energy distributions of both electron and muon events, the 
SuperK experiment convincingly established that the muon-neutrino flux has a deficit 
compared with the calculated flux that increased with zenith angle (or equivalently 
the path distance), while the electron- neutrino flux agreed with no-oscillation ex- 
^For a historical perspective and discussions of early work on neutrinos, see Refs. 
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pectations. The accumulation of high statistics data by SuperK eventually excluded 
interpretations other than z/^ to z/^ oscillations with nearly maximal mixing at a mass- 
squared- difference scale ~ 2.0 x 10"'^ eV^. The energy and angular resolution of 
the SuperK experiment are too coarse to allow the first minimum in z/^ —>■ to be 
resolved. That is the goal in the ongoing K2K jTHj and the forthcoming MINOS J7] 
and CNGS JHl UH] accelerator experiments. 

The evidence for solar neutrino oscillations continued to build as experiments with 
different detectors and energy sensitivity all found deficits of 1/3 to 1/2 in rates com- 
pared to the SSM. The ''^Ga radiochemical experiments of SAGE 20j, GALLEX [21] 
and GNO [221 found deficits in the flux of pp neutrinos, which is the dominant prod- 
uct of the pp reaction chain that powers the Sun^. The SuperK water Cherenkov 
detector [24] accurately measured the electron energy spectrum from high-energy so- 
lar neutrinos (>5 MeV) originating from decays in the Sun and found it to be 
flat [23 I2n] with respect to the SSM prediction. A crucial theoretical aspect in the 
oscillations of solar neutrinos is coherent forward u^e —>■ u^e scattering, first discussed 
by Wolfenstein j^, which affects electron neutrinos as they propagate through the 
dense solar core. Matter effects can produce large changes in the oscillation amplitude 
and wavelength compared to vacuum oscillations, as first shown by Barger, Whisnant, 
Pakvasa and Phillips j2H]- A resonance enhancement^ can be realized in matter for 
one sign of 5m^. Because of the prevailing prejudice that neutrino mixing would 
be small, there was a strong theoretical bias in favor of a resonant solar solution, 
which was the original solution to the solar neutrino problem proposed by Mikheyev 
and Smirnov (the so-called Mikheyev-Smirnov- Wolfenstein or MSW solution) [21] . In 
addition to this small mixing angle solution (known as SMA), other solutions with 
a large vacuum mixing angle were later identified that could account for the solar 

■^The pp neutrino flux is now determined experimentally to ±2% and is in agreement with the 
SSM predictions to 1% (the theoretical uncertainty is also 1%) [22] i thus confirming the essential 
ingredients of the SSM. 

•^BWPP discovered this enhancement in a medium of constant density by varying the neutrino 
energy; the matter effect depends on the product of the neutrino energy and the electron density. 
Mikheyev and Smirnov applied the enhancement at a given neutrino energy to the propagation of 
solar neutrinos through the varying electron density in the Sun. 
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neutrino flux suppression [HUj . 

The other solutions were named LMA (large mixing angle), LOW (low Sm^, 
low probability) QVO (quasi- vacuum oscillations) and VO (vacuum os- 

cillations) jHSl- These islands in the (5m^,tan^^) in the solar neutrino oscillation 
parameter space are illustrated in Fig. ^ The fiat energy spectrum relative to the 
SSM and the absence of a significant day/night difference (that can be caused by 
Earth- matter effects), measured by SuperK (see Fig. |21), led to a strong preference 
for the LMA solution. Subsequently, the Sudbury Neutrino Observatory (SNO) [33] 
measured the neutrino neutral currents as well as the charged currents, and confirmed 
the oscillations independent of the SSM normalization of the flux. The recent SNO 
salt phase data ^| in conjunction with other solar neutrino data selects the LMA 
solution uniquely at more than the 3a level. The mass-squared difference indicated 
by the solar neutrino data is ~ 6 x 10^^ eV^ and the mixing is large but not maximal, 
tan^ 9 ~ 0.4; see Section 4. 




Iff^ 10'^ 10'^ 10'^ 1 10. 



tan^G 

Figure 1: 90%, 95%, 99% and 3a C. L. allowed regions in the (5m^, tan^ 9) oscillation 
parameter space, before any SNO data. From Ref. |34j . 



6 



^0.75 
n 7 
|o.65 

CO 

Q 0.6 
0.55 
0.5 
0.45 



10 


-10 
-20 
-30 
-40 
-50 



LMA best-fit Spectrum 



— I — I — I — I — I — I — I — I — I — — i-^-i — 1,1,1 — I — I — I — I — I 



toy. 



±1o 



LMA best-fit D/N Asymmetry, 



6 8 10 12 14 16 18 20 

E in MeV 



Figure 2: The energy spectrum and the day-night asymmetry (day rate — night 
rate divided by the average rate) as measured by SuperK. The sohd hues are the 
predictions for 6ml = 6.3 x 10^^ eV^ and tan^ 6'^ = 0.55. From Ref. |26|. 



The definitive confirmation of the LMA solar solution has come from the Kam- 
LAND experiment [HZj- If CPT invariance is assumed, the probabilities of z/g — > Ue 
and Ue — i> z/g oscillations should be equal at the same values of L/E. At the average 
distance L ~ 180 km of the reactors from the KamLAND detector and the typical 
energies of a few MeV of the reactor z/g, the experiment has near optimal sensitivity to 
the (5m^ value of the LMA solar solution. The first year of data from KamLAND [SB] 
shows the rate suppression of P(z/e —>■ i^e) expected from the solar LMA solution, 
vindicating the oscillation interpretation of the solar neutrino problem. The solar 
are oscillating to a combination of z/^ and Ur. 

With three neutrinos (z/g, z/^, z/^) there are two distinct (5m^ that can account for 
the atmospheric and solar neutrino oscillations. However, the LSND accelerator ex- 
periment [Sn] found evidence at 3.3cr significance for z/^ —>■ oscillations at a higher 
5rn? ~ 0.2 to 1 eV^ with small mixing, sin^26' ~ 0.003 to 0.04. Evidence for os- 
cillations of — > Ve was also observed, although at lesser significance [IHl- Such 
oscillations are constrained but not excluded by the KARMEN experiment jJIJ- A 
possible way to explain the LSND effect is to invoke oscillations involving a fourth, 
sterile neutrino with no Standard Model interactions 02] • Two possible scenarios are 
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considered, with 2 + 2 and 3 + 1 spectra of (5m^ jlHl HH ESI HE] , where the large 
mass-squared difference of the LSND oscillations connects the upper and lower levels. 
However, there is a tension in the global fits between the solar and atmospheric data 
on one hand and the reactor and accelerator data on the other, with the conclusion 
that both schemes are highly disfavored |17]. The MiniBooNE j3H] experiment, now 
in progress at Fermilab, is designed to test the u^. oscillation hypothesis for the 
LSND effect. The source of the antineutrino beam at LSND is muon decays, while at 
MiniBooNE the source is pion decays. An alternative nonoscillation interpretation of 
the LSND effect as a nonstandard muon decay ^1 will not be tested by MiniBooNE. 
Since the LSND effect is fragile, this review emphasizes three neutrino oscillation 
phenomena. 

The phenomenology of three neutrino mixing, like that of three quarks in- 
volves three mixing angles {623 = 9a, 612 = 9s, 613 = 9^), two mass-squared differences 
(^mgi = 6ml, ^"^21 = ^^l)^ and one CP- violating phase (6) [3 EH E21 ESI EH ■ In 
order that 6 be measurable, both scales must contribute to the oscillation |52j . 
If neutrinos are Majorana ^j, two further CP- violating phases (02, 03 ) enter in 
the calculation of neutrinoless double-beta decay [56J but not oscillations [S^. We 
already have approximate knowledge of 9 a, 9s, 6m1 (but not the sign of Sml^) and 
(its sign is known from solar matter effects). The major challenge before us now is 
the measurement of 9x, for which we have only an upper bound, sm'^29x ^ 0.2 (for 
6ml = 2.0 X 10"=^ eV2) at the 95% C. L. from the CHOOZ [58; and Palo Verde ^ 
reactor experiments. 

The amplitude of the oscillations z/^ z/g and u^. —>■ Uj. are governed by the size of 
9^. Long-baseline and reactor experiments are planned with improved sensitivity to 
sin^2^3; innj; see Section 10. The CP- violating phase enters oscillations via a factor 
sin6'a;e~*'^ and thus it is essential to establish a nonzero value of 9^ in order to pursue 
the measurement of 6. The size of CP violation in long-baseline experiments also 
depends on the value of 6ml. Oscillation parameter ambiguities exist [HI] that must 
be resolved by the experiments. 

The anticipated steps in the long-baseline program are off-axis beams [UmiH^ IUS], 
superbeams jHH ES] , detectors with larger fiducial volumes and sophistication [^01 
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inniinZI, and neutrmo factories jHElinni- The ultimate sensitivities will be derived from 
neutrino factories, where the neutrino beams are obtained from the decays of muons 
that are stored in a ring with straight sections. 

Although oscillations have allowed us to establish that neutrinos have mass, they 
do not probe the absolute neutrino mass scale. In particle and nuclear physics, the 
only avenues for this purpose are tritium beta decay and neutrinoless double-beta 
decay, and the latter only if neutrinos are Majorana particles. Experiments are be- 
ginning to probe the interesting region of m,y. Another route to the absolute mass is 
the power spectrum of galaxies, which gets modified on small scales when the sum of 
neutrino masses is nonzero The WMAP analysis [71] of the cosmic microwave 
background (CMB) and large-scale structure data have already given an upper limit 
on about 1 eV. 

The field of neutrino physics is progressing at a rapid rate. The purpose of this 
review is to summarize the present status of the field and to discuss ways that progress 
can be made in experimentally answering the outstanding questions. For other recent 
reviews see Refs. [721 ESI IZH ESI EO] • We have provided an extensive bibliography, 
but due to the large number of papers on the subject, it may not be comprehensive. 



2 Neutrino counting: Z-decays, CMB and BBN 

Studies of e^e~ annihilation at the Z-resonance pole at the Large Electron Positron 
collider have determined the invisible width of the Z boson. The experimental value 
Ny = 2.984 ± 0.008 is close to the number expected from 3 active light neutrinos, 
though the value is 2cr low [54J. 

The cosmic microwave background (CMB) anisotropies and Big Bang Nucleosyn- 
thesis (BBN) probe the effective number of neutrinos {N^, = 3 + AA^,^) that were 
present in the early universe. The extra relativistic energy density due to sterile neu- 
trinos, or other possible light particles, is normalized to that of an equivalent neutrino 
flavor as [77] 

= AN.p, = l^Kp^ , (1) 

o 

where p^ is the energy density in photons. Sterile neutrinos would contribute to AA'y, 
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but so could other new physics sources. 

The precise WMAP measurements jTJ of the CMB have been analyzed to con- 
strain AA^,^ [7H1 ESI • f^^f • EH] , a flat universe with a cosmological constant A as 
dark energy is assumed. The parameters varied are the reduced Hubble constant 
the baryon density ujb = ^sh'^, the total matter density Um = ^hih^ (comprised 
of baryons and cold dark matter), the optical depth r, the spectral index Ug and 
amplitude As of the primordial power spectrum (P = As{k / k^)"""'^) and AA^^. The 
HST measurement of h = 0.72 ± 0.08 'HDI is imposed as a top-hat distribution. This 
strong h prior helps to break the degeneracy between um and AA^^. The universe 
is also required to be older than 11 Gyr, as inferred from globular clusters j81j . 
The effects of on the CMB are illustrated in Fig. El for a fixed normalization 
of the power spectrum. The resulting AA^i, constraints are shown in Fig. |31 where 
r/io = lO^^UB/n^ = 21Aujb. The best-fit is AA^^ = -0.25, but at 2ct, AA^^ < 5.3 is 
allowed; see Table [T] [78j . 




Multipole 1 

Figure 3: The power spectrum for the best-fit {N^ = 2.75) to the WMAP data [7T] 
is the solid line. With all other parameters and the overall normalization of the 
primordial spectrum fixed, the spectra for A^^ = 1, A^j, = 5 and = 7 are the 
dotted, dot-dashed and dashed lines, respectively. The data points represent the 
binned TT power spectrum from WMAP. From Ref. |78j . 
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Figure 4: The la and 2cr contours in the rjiQ — AN^, plane from an analysis [TH] of 
WMAP data. The solid (dotted) lines correspond to to > H (12) Gyr. The cross 
marks the best-fit at uj-q = 0.023 and AA^,^ = —0.25. 

Table 1: The 2a ranges (for 1 degree of freedom) of N^, and r]io from analyses [TH] 
of WMAP data, deuterium and helium abundances and their combinations. The 
WMAP analysis involves the assumption of a flat universe, along with the strong 
HST prior on h and the age constraint to > 11 Gyr. For BBN the adopted primordial 
abundances are: ?/d = 10^(D/H)= 2.6 ± 0.4 [H2], Y = 0.238 ± 0.005 [H3]- 





Ny {2a range) 


?7io {2a range) 


WMAP 


0.9 - 8.3 


5.58 - 7.26 


VB + Y 


1.7-3.0 


4.84 - 7.11 


WMAP + yD+Y 


1.7-3.0 


5.53 - 6.76 



BBN is a much better probe of AA^i^ than the CMB. The prediction of the pri- 
mordial abundance of ^He depends sensitively on the early expansion rate, while the 
prediction of the D abundance is most sensitive to the baryon density rjiQ [H^. The 
best-fit value of AA^^ from BBN is AA^^, = -0.7 but at 2a, A"^ = 3 is allowed. The D- 
inferred baryon density is in excellent agreement with the baryon density determined 
from the CMB and Large Scale Structure. 
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A combined analysis of the BBN and the WMAP data yields the allowed regions 
in Fig. El As with the BBN analysis above, the best-fit value is /S.Ny = —0.7 and 
A^;^ = 3 is allowed at 2a. 



<] 




Figure 5: The la, 2a and 3a contours in the rjiQ — AN^, plane from a combination of 
WMAP data and the adopted D and "^He abundances. From Ref. [78] . 



The BBN analysis is consistent with 3 neutrinos, but gives no support to the 
possible existence of extra neutrinos. Even one extra, fully thermalized neutrino 
{ANi, = 1) is strongly disfavored [75|. We discuss the implications of this result for 
the sterile neutrino interpretation of the LSND experiment in Section ITTl 



3 Neutrino mixing and oscillations 
3.1 Vacuum oscillations 

The dramatic increase in our knowledge of neutrino properties has come from obser- 
vational evidence of neutrino oscillations. These neutrino flavor changes require that 
the neutrino flavor states, z/^ are not the same as the neutrino mass eigenstates, z/j. 
The eigenstates are related by a unitary matrix V ff. , 

^a = Y.^>^- (2) 
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V is often denoted as Vmns, where MNS represents the authors of Ref. [7j^. For 
3 neutrinos, the mixing matrix V is specified by three rotation angles 6a, 0^, Og 
< 6i < 7r/2) and three CP-violating phases 5, 02 and 03 (0 < 5, 0i < 2tx). V can 
be conveniently written as the matrix product 



V 



1 







Ca 





Sa 
Ca. 








SxC 





1 










Cx 



Cs 





Ss 
Cs 








1 



1 

e"V2 













(3) 

where q denotes cos 6*4 and Si denotes sin 6*1. The angle 9a, customarily denoted as 6*23, 
governs the oscillations of atmospheric neutrinos, the angle 6s (^12) describes solar 
neutrino oscillations, and the angle 6x {613) is an unknown angle that is bounded 
by reactor neutrino experiments at short distances (L ~ 1 km). The oscillation 
probabilities are independent of the Majorana phases 02 and 03. Vacuum neutrino 
oscillations are given by 

2 



(4) 



where the are the neutrino eigenmasses. The oscillation probabilities depend 
only on differences of squared neutrino masses. The oscillation arguments for the 
atmospheric and solar phenomena are 



brr?aL 

Aa = — 



respectively, where 



1 ) 



rrin 



mi 



Srria = m\ — m 
The vacuum oscillation probabilities are 

P(z/e^i/e) = 1 - sin^ sin^ Aa - (c^ sin^ 2^, + sin^ 2^x) sin^ A, 



(5) 



(6) 



+.^in^2^.(-sin2A.sin2A. + 2sin^A.sin^A. 

„2 „-„2 



(7) 



P{ve^i^ti) = sin^ 26'x sin^ Aa + 4J(sin2A^ sin^ Aa - sin2Aa sin^ As) 
^Sometimes it is denoted as Vpmns or Vmnsp to acknowledge the contributions of Pontecorvo. 
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-{slsl sin^ 26.^ - AK) [- sin 2 A, sin 2A„ + 2 sin^ A, sin^ A„] 
+ [cl{cl - slsl) sin^ 26, + sis] sin^ 26, - SKsl] sin' A, , (8) 
P{i^^^iy^) = l-{4sm'26a + slsm'26,)sm'Aa 

+[cl{cl - slsl) sin^ 26a + s^sl sin^ 26, - 8Ksl] 

X sin 2 As sin 2Aa + 2 sin^ A^ sin^ A^] 
-[sin^ 26sicl - slsiy + si sin^ 20„(1 - cj sin^ 26 

+2sx sin 26s cos 26s sin 26'a cos 26aCs — WKs^sl 

+ sin^ 26acl{cl - slsl) + slsl sin^ 26,] sin^ A, (9) 
Pi^ii-^i^r) = 4sin^26'asin^ Aa + 4J(sin2AsSin^ Aa -sin2AaSin^ As) 
-[cl sin' 26a{cl - slsl) + AK cos 26 a] 

X sin 2 As sin 2 A^ + 2 sin^ A^ sin^ A J 
+ [sin2 26aicl - slsl)' + si sin' 20^(1 - sin^ 26acl) + AKcos26a 

+s, sin 26s cos 26's sin 26a cos 26a{l + sl)cs] sin' As , (10) 

where the CP-violating quantity J |HS1 is 

J = -Ca; sin 26*3; sin26'a sin 26'sS5 , (11) 
8 

and for convenience we have defined 

K = -c, sin 26, sin 26a sin 26 sC^ . (12) 
8 

Oscillation probabilities for other neutrino channels may be obtained by probability 
conservation, z.e., ^^^Piya ^ ^/s) = J^fs^i^a ^p) = 1- Probabilities for an- 
tineutrino channels are obtained by replacing 5 by —5 in the corresponding neutrino 
formulae. Also, since CPT is conserved for ordinary neutrino oscillations, the an- 
tineutrino probabilities are given by P{pa i^p) = P{^i3 ^a), and P(z/a I'/s) is 
obtained from P(z//3 — > u^) by replacing 6 by —6. 

For oscillations of atmospheric and long-baseline neutrinos, the oscillation argu- 
ment Aa is dominant, and the vacuum oscillation probabilities in the leading oscilla- 
tion approximation (where only the 6ml oscillations are appreciable) are 

P(z/, -.z/,) ~ l-sin2 2^,sin2A„ (13) 
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~ si sin^ 26.^ sin^ 

~ 1 - (c^ sin^ 29 a + si sin^ 29^) sin^ 
~ 4 sin^ sin^ A„ . 



(14) 
(15) 
(16) 



For vacuum oscillations of solar neutrinos, where |Aa| ^ 1, the terms involving 
Aq approach their average values over a complete cycle and 

P(z/e z/e) ~ 1 - ^ sin^ 29.^ - c\ sin^ 29, sin^ A, . (17) 

Note that in the limit 9^ — »■ 0, z/^ — > oscillations of atmospheric neutrinos and 
— > oscillations of solar neutrinos completely decouple, i.e., they are determined 
by independent parameters, and each has the form of two-neutrino oscillations. 



3.2 Matter effects on oscillations 

The scattering of v^. on electrons in matter can modify the vacuum oscillation proba- 
bilities j2Zl I2HI EHl • For the two-neutrino case, with mixing angle 9 and mass-squared 
difference brv? , the oscillation probability amplitude sin^ 29 m in matter is 

9 . sin2 2^ , , 

sin' 29^ = — -2 — , 18 

(^-cos2^) +sin2 29 

where 

A = 2V2 GpN.E^ = 1.54 X 10"^ eV^ FeP(g/cm^) E^(GeV) , (19) 

and Ne is the electron number density, which is the product of the electron fraction 
Ye and matter density p. The oscillation amplitude in matter is enhanced if > 0, 
and a resonance occurs {i.e., the amplitude reaches its maximal value of unity) at 
the critical density = 5m? cos 29 / {2\/2Gf E^). For antineutrinos, A — > —A in 
Eq. p8|) . and the oscillation amplitude in matter is enhanced if 5ir? < 0. For A^^e 
much larger than the critical density, the oscillation amplitude in matter is strongly 
suppressed for both neutrinos and antineutrinos. The effective value of 5rr?, and 
hence the oscillation wavelength, is also changed in matter: 




5mi^ = 5m\\\^^- cos 29] +sin^2^. (20) 
15 







COS 9m — sin 9m 










sin 9m COS 9m 







The angle 9m represents the mixing between the flavor eigenstates and the 
instantaneous eigenstates in matter Vim'- 



(21) 



If the electron density is iVg when the neutrino is created, the initial value of 6*^ is 
given by 

„ - cos 26' , ^ 

cos 291. = , (22) 

■ AO 

■ 5m'^ 



COS 2^)2 + sin2 29 

where A" = 2^2 E^. A neutrino originally created as a z/g can be expressed in 

terms of the lower and upper eigenstates as z/g = cos 9m i^im + sin 9m i^2m- 



3.3 Solar neutrino oscillations 

For matter of varying density, such as for neutrinos propagating through the Sun, the 
instantaneous eigenstates change as the neutrinos propagate. Once a solar neutrino 
reaches the vacuum, the lower eigenstate is ui = cos 9Vf. — sin 9 and the upper 
eigenstate is V2 = sin 6' z/g + cos 9 v^. 

A z/g created far above the resonance density is predominantly in the upper eigen- 
state. Then if the neutrino propagation is adiabatic, the neutrino will remain in the 
upper eigenstate, and if 9 is small it will be predominantly z/^ once it reaches the vac- 
uum [SIIIHH] (see Fig. ini). For nonadiabatic propagation, if the probability of jumping 
from one eigenstate to another is P^, then averaging over the oscillations 



0) = - [1 + (1 - 2P,) cos 29 cos 2^] 



where jHOl 



exp( - I 7F) - exp( 



2 ' sir?6 



(23) 



(24) 



l-exp(- 27-?,. 

is the transition probability with P = 1 — tan^^ for the exponentially varying matter 
density in the Sun, and [88, 
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4:V2GfEI \dNe/dLl 
16 



(25) 




vacuum 



core 



Figure 6: The level-crossing diagram for solar neutrinos. 

measures the degree of adiabaticity of the transition. In Eq. ()25j) . \dNe/dL\c is the 
density gradient at the critical density. This process is analogous to level crossings in 
atoms [HI] . For an electron neutrino that is created well above the resonance density 
[9^ ~ 7r/2) and which undergoes a perfectly adiabatic transition (7 — > 00, — > 0), 
the oscillation probability is P(z/e Ue) = sin^ 9. Thus a very large depletion of solar 
z/g's is possible even for small vacuum mixing angles. This is known as the MSW 
effect, and was first studied numerically in Ref. ^]. In the extreme nonadiabatic limit 
(7 — > 0) [92J, Px — cos^ 9 and the oscillation probability approaches 1 — | sin^ 29, the 
expected value for two-neutrino vacuum oscillations averaged over the oscillations. 



The range over which P{i'e ^e) < h 



IS 



< El, < 6m'' sm 29 J — — . (26) 



2V2GfN^ \ 8V2\n2GF\dNe/dL\ 

For neutrino energies below this range, the initial density is below the critical density 
and the neutrino starts with a large fraction in the lower eigenstate. For neutrino 
energies above this range, the transition becomes very nonadiabatic and the neutrino 
has a high probability of hopping from the upper eigenstate to the lower eigenstate. 
In either large component of the neutrino ends up in the lower eigenstate, 

in which case the survival probability is greater than |. For = lOOA^^/cm^ 
(the approximate number density at the center of the Sun), neutrinos with energies 
in the range 2 MeV < Eu < 20 GeV will have survival probabilities smaller than |, 
assuming the best-fit oscillation parameters of the LMA solution {Smi = 7 x 10^^ eV^ 
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and sm'^2es = 0.83). 

Exact formulae that include cases when the neutrino is created near resonance 
are presented in Ref. jHS]- More discussions of exact formulae for the transition 
probability are given in Ref. [M]. A semi-classical treatment for an arbitrary density 
profile is given in Ref. j^Sl- Formulae for the MSW effect in a three-neutrino model 
are presented in Ref. [HB] . 



3.4 Long-baseline oscillations through the Earth 

Oscillations of long-baseline neutrinos are affected by electrons in the Earth if the 
path length is an appreciable fraction of the Earth's diameter. The full propagation 
equations for three neutrinos in matter are 



2K^ 



(27) 



A constant density approximation often provides a good representation of the neu- 
trino propagation over long baselines through the Earth. Since -C and Qx 
is small, the probabilities can be expanded to second order in terms of the small pa- 
rameters Qx and |5m^/5m^| jH7|inH]- The following useful approximations for > 
are obtained [HT] : 



P + 2xyfg{cos 5 cos A — sin 5 sin A) + y'^g'^ , (28) 
x'^P + 2xyfg(cos S cos A + sin 5 sin A) + y^g"^ , (29) 



where 



and 



X 

y 

fj 

9 



A= |A„ 



sin 9a sin 29^ , 
a cos 9a sin 29s 



sm 



lTi)A /(iTi), 



sin(iA)/i, 



A=\A/Sml\ 



a = \5ml/5ma\ 



(30) 
(31) 
(32) 
(33) 

(34) 
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For Sm^ 



< 0, the corresponding formulae are 




^2 j2 _ 2xyfg(cos S cos A + sin 5 sin A) + y^g'^ , 



(35) 




^2 j2 _ 2xyfg{cos S cos A — sin 5 sin A) + y'^g^ . 



(36) 



Oscillation probabilities for an initial and final z/^ can be found by changing the 
sign of the sin 5 term in Eqs. (PHl EHl ESI), and (jHSl)- These expansions are nearly- 
exact for distances less than 4000 km when Ei, > 0.5 GeV [61J. For more accurate 
results at longer distances, Eq. pTjl may be integrated numerically over the density 
profile of the neutrino path. 

Other approximate solutions can also be useful in certain situations. At relatively 
short distances where the matter effect is not as large, an expansion can be made 
in a and A/6m1 for the constant density solution [221 • Relationships between the 
vacuum and matter oscillation parameters for three-neutrino oscillations are given 
in Ref. |100j . Exact results for the three- neutrino case with constant density are 
given in Refs. [2H1 llOlj . Several properties of the general three neutrino solution with 
a nonconstant density profile are discussed in Ref. |l()2j . Consequences of random 
density fluctuations are discussed in Ref. jl()3j : they are not expected to play an 
important role in most situations. 

The evolution in Eq. (j2Z|) is modified if a sterile neutrino Ug, is involved |104j : 



where Nn is the neutron number density. In two-neutrino oscillations between 
and a sterile neutrino, the electron number density N^. is changed to an effective 



in oscillations of solar neutrinos. Also, in two-neutrino oscillations between a or Uj. 
and a sterile neutrino, iVg// = — |A^„; consequently there can be substantial matter 
effects in z/^ —>■ Ug oscillations of atmospheric neutrinos. 



4 The solution to the solar neutrino problem 

Decades of study of neutrinos from the Sun P 1201 IHl 122 123 123 ESI EE] have con- 
vincingly established that neutrino oscillations are the cause of the deficits of 1/3 to 




(37) 



number density N^ff = Nf, — ^ 



^Nn, which changes the critical density for a resonance 
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1/2 in the measured electron- neutrino flux relative to the Standard Solar Model ex- 
pectations ^T]. The water Cherenkov experiments of SuperK and SNO measure the 
high energy neutrinos {E > 5 MeV) from the *B chain, the Chlorine experiment also 
detects the intermediate energy neutrinos from ^Be and pep, and the GALLEX, GNO 
and SAGE experiments have dominant contributions from the pp neutrinos [TOj; see 
Fig. [7| . Until recently, the interpretation of the deficits depended on comparisons 
with SSM predictions of the flux. With the SNO experiment, which directly measures 
the total active neutrino flux via neutral currents, the evidence for flavor conversion 
becomes robust. 




Neutrino energy (MeV) 

Figure 7: The neutrino flux predictions of the Standard Solar Model ^T\. From 
Ref. ISH. 

The SNO experiment utilizes a heavy water target and measures the following 
processes jHSl EH] '■ 

Charged- Current (CC): i'e + d^e~+p + p (38) 
Neutral- Current (NC): + d ^ u,^ + n + p (39) 

Elastic-Scattering (ES): z/^. -|- e~ — + e~ (40) 

The CC/NC ratio establishes the oscillations of z/g to and u^- flavors, 

CC ^ flux(z/e) , . 

NC flux(l/e + ly^^ + lyr)' 
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Only z/e are produced in the Sun; the z/^ and z/,- fluxes are a consequence of oscillations. 
The charged- current signal was found to be suppressed by 7.6a from the neutral- 
current signal; see Fig. |H1 



Neutral Elastic Charged 

Current (NC) Scattering (ES) Current (CC) 




Figure 8: Evidence for neutrino flavor change seen by SNO. The open (filled) circles 
represent the 2003 SNO flux results, relative to the SSM, under the assumption of an 
undistorted (unconstrained) ®B neutrino energy spectrum. 

The day and night energy spectra of charged- current events are potentially sen- 
sitive to matter effects on oscillations that occur when the neutrinos travel through 
the Earth |105j . though a significant day-night asymmetry is yet to be established. 
SuperK (SNO) measured the day minus nig ht Ue rate to be -2.1% ± 2.0%^};^% 
(-7.0% ± 4:.9%t\i%) of the average rate |2SII2S] (EH]); see Fig. M From a global fit 
to neutrino data, regions of the solar oscillation parameters have been determined, as 
shown in Fig. . The Large Mixing Angle (LMA) solution is preferred at more 

than 3cr. The best fit to the solar data is (5m^ = 6.5 x 10"^ eV^ and tan^ 6s = 0.40 jHEj- 
The survival probability versus neutrino energy for LMA parameters is shown in 
Fig.im 

The KamLAND experiment jSZI measures the electron antineutrino flux at the 
Kamiokande site from surrounding reactors. The dominant reactor is at L = 160 km 
and the average distance from the sources is L ~ 180 km. The measured reaction is 
^'e + P ~* + n. If CPT invariance holds, which is expected in quantum field the- 
ory, then P{ue ^e) = P{ve i^e)- Therefore, for the LMA solar solution, reactor 
antineutrinos should also disappear due to oscillations. For any other solar oscilla- 
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5 10 15 20 25 

5m^ (10'W) 



Figure 9: The SuperK day-night asymmetry as a function of Srnl, with the la band 
shaded. The sohd hne is the prediction for tan^^^^ = 0.55, and the cross-hatched 
bands are the ranges of Sml allowed by KamLAND. From Ref. ^B] . 
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Figure 10: 90%, 95%, 99% and 3cr C. L. allowed regions from a fit to the Homestake, 
GALLEX+GNO and SAGE rates, and the SuperK and SNO spectra (with NC sensi- 
tivity enhanced by salt). The ellipse is a projection of the 3a region from three years 
of KamLAND data assuming the best-fit LMA parameters {Sm^ = 6.5 x 10~^ eV^, 
tan^^, = 0.40). Adapted from Ref. ^ and Ref. [TUHj . 

tion solution, no disappearance would be observed at KamLAND. The KamLAND 
expectation |106| I108j for 3 years data, assuming the LMA oscillation parameters, is 
shown in Fig. by the narrow ellipse superimposed on the LMA region from solar 
data. With sufficient data, the KamLAND experiment should "see" the oscillations 
in the positron energy spectra, as illustrated in Fig. [T21 |106j . 
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Figure 11: The flux- weighted survival probability for an LMA solution. Adapted 
from Ref. ITOTj . 



300 



100 



No Oscillations 

/"^ 5m| = 7 X 


10"^ 




/ \ 5m| = 1 X 


lO-' 


eV^ 


/ \ sin^aSs = 


0.75 




/ ' ' \ 
It '» ' ^ \ 

/ ■' ■/' \ 

/ / A 'A 

/ ■' /' * \ 
/ I ' \ * \ 






\ \ * \ 
/ 'i ■ » \ 

/ *■ \ \ 

- //'' \ 

/ ' »' \ 

/'; *v \ 






/'■ * ^ \ 
/(* * * \ 

/f ^ \ 

/' A 






- /'' \v 
/ 







) 3 4 6 
Total e* Energy (MeV) 




8 



Figure 12: KamLAND is expected to detect the spectral distortion resulting from 
oscillations. Adapted from Ref |lU6j . 
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The first KamLAND results [38] are based on 145 days of operation. Tlie data 
give spectacular confirmation of the solar oscillation analysis predictions; see Fig. IT^ 
The fractional number of events (A^(observed) — A^(bkg))/A^(expected) = 0.611 ± 
0.085(stat) ± 0.041 (syst) excludes no oscillations at 99.95% C. L. and eliminates all 
solar solutions but LMA PB] • A combined analysis of solar and KamLAND data select 
the LMA solution uniquely at the 4 — 5o" C. L |110j . Maximal mixing is excluded at 
the 5.4a C. L. and 5mf < 10"^ eV^ at greater than the 99% C. L. [3^1; see Fig. d 
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Figure 13: The ratio of the measured to expected flux from reactor experiments. 
The shaded region encompasses the fluxes predicted for oscillation parameters in the 
95% C. L. region from an analysis of pre-SNO salt phase data |109j . The dotted curve 
corresponds to Smj = 5.5 x 10"^ eV^ and tan^ 6s = 0.42. From Ref. jSH]- 



KamLAND has not only eliminated all oscillation solutions other than LMA, but 
has relegated nonoscillation solutions to the solar neutrino problem to be at most 
subleading effects. Non-standard neutrino interactions (NSNl) [1111, which lead to 
energy-independent conversion probabilities, were consistent with the flat energy spec- 
tra seen by SuperK and SNO jll2j . With KamLAND data, NSNI are generically 
rejected as the leading cause of z/g-disappearance at about the 3cr C. L. Also, the 
resonant and nonresonant spin-flavor precession solutions |113j are allowed only at 
the 99.86% and 99.88% C. L., respectively |114j . Yet another excluded alternative 
invokes the violation of the equivalence principle to induce oscillations even for mass- 
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Figure 14: 90%, 95%, 99% and 3cr C. L. allowed regions from a combined fit to 
KamLAND and solar neutrino data. The best-fit point is at 5mf = 7.1 x 10"^ eV^ 
and tan^^, = 0.41. From Ref. |SE]. 

less neutrinos |115j . These three solutions fail because the KamLAND baseline is too 
short for any significant disappearance to occur. 

The continuation of the KamLAND reactor experiment will provide a measure- 
ment of 6m1 that will be precise to about 10%. Data from solar neutrinos require that 
the sign of 6m1 is positive and that the mixing angle 6s is nonmaximal. A more precise 
determination of 6s is important to test models of neutrino mass. Future SNO data 
(from the ^He proportional counter phase) should reduce the presently allowed range 
of 6s by measuring the CC/NC ratio and the day-night asymmetry more precisely. As 
far as the parameters responsible for the mixing of solar neutrinos are concerned, the 
goal of any new solar or reactor experiment should be to pin down 6s more accurately 
than a combination of future SNO and KamLAND data. Borexino |116j or any other 
^Be solar neutrino experiment will not improve on the accuracy with which three 
years of KamLAND data and solar data will determine 6s [23^ . A future pp solar neu- 
trino experiment with better than 3% precision can lead to significant improvement. 
Proposals for measuring pp neutrinos include LENS |117j . MOON |118j . SIREN |119j . 
XMASS [1201, CLEAN [Ml, HERON [1221 and GENIUS [123. Another avenue for 
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an improved 6^ measurement is a lithium-based radiochemical detector that detects 
neutrinos from the CNO cycle |124j . A reactor neutrino experiment with baseline 
such that the measured survival probability is a minimum can lead to a more precise 
measurement of Og [125J. Since the reactor neutrino spectrum is accurately known 
and KamLAND will determine precisely, such an experiment is conceivable. For 
Sm^ = 7 X 10~^ eV^, the required baseline is 70 km. 

5 Atmospheric neutrinos 

The first compelling evidence for neutrino oscillations came from the measurement of 
atmospheric neutrinos. Interactions of cosmic rays with the atmosphere produce pions 
and kaons that decay to muon neutrinos, electron neutrinos, and their antineutrinos: 

7r+, z/^./i+ z^^e+z/eZ/^ , (42) 

7r~,K'' v^iT z/^e~z/ei/^ . (43) 

On average there are twice as many muon neutrinos as electron neutrinos at energies 
of about 1 GeV, although the electron neutrinos tend to be at somewhat lower energies 
since they are produced only in a secondary decay. The atmospheric neutrino flux is 
well understood: the normalizations are known to 20% or better and ratios of fluxes 
are known to 5% |126j . The flux falls off rapidly with neutrino energy for > 1 GeV. 

Neutrinos observed at different zenith angles have path distances that vary from 
L ~ 10 — 30 km for downward neutrinos to L ~ 10^ km for upward neutrinos, as 
illustrated in Fig. ^| The ratio of observed to expected neutrino events provides a 
sensitive measure of neutrino oscillations, especially since different values of neutrino 
baselines and energies can be studied. Initial evidence for atmospheric neutrino oscil- 
lations was an overall depletion of muon neutrinos compared to the theoretical 
expectation. The SuperK experiment [T3j has studied CC events in four categories: 
fully contained {E^, ~ 1 GeV), partially contained {E^, ~ 10 GeV), upward-going 
stopped {El, ~ 10 GeV) and through-going {E^, ~ 100 GeV). The contained events 
have the highest statistics and give the more precise measurement, but all of the data 
samples are fully consistent with the same oscillation parameters. Zenith angle dis- 
tributions for the e-like and /x-like contained events are shown in Fig. along with 
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Figure 15: A schematic view (not to scale) of the different zenith angles of atmospheric 
neutrinos and distances they travel before detection. 

the no oscillation expectation and the best fit assuming oscillations. Atmospheric 
neutrinos are primarily sensitive to the leading oscillation which involves 9a and Ox- 
The angle Ox is constrained to be smaller than 13° (for 5m?^ = 2.0 x 10^^ eV^) at the 
95% C. L. Assuming 6'^^ = in Eqs. and (fTUI) {i.e., atmospheric muon neutrinos 
oscillate exclusively to tau neutrinos), the latest (preliminary) analysis by the SuperK 
collaboration of their data yields best-fit values sin^ 20a = 1-00 (maximal mixing) and 
= 2.0 X 10~^eV^ |127j . The 90% C. L. ranges for these oscillation parameters 
are sin^ 20„ > 0.90 and Sml ~ (1.3 - 3.0) x lO'^eV^. Both the Soudan-2 fTI^ 
and MACRO |129j experiments have also measured atmospheric neutrinos and find 
allowed regions consistent with the SuperK result (see Fig. I17|l. 
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Figure 16: Zenith angle distributions for e-like and /x-like contained atmospheric 
neutrino events in SuperK [127^ ; cos = 1 corresponds to downward events with L ~ 
15 km and cosG = —1 corresponds to upward events with L ~ 13000 km. The hues 
show the best fits with and without oscillations; the best-fit is 5m?^ = 2.0 x 10~^ eV^ 
and sin^ 2^„ = 1.00. 
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Figure 17: 90% C. L. allowed regions for z/^ oscillations of atmospheric neutrinos 
for Kamiokande, SuperK, Soudan-2 and MACRO. From Ref. |130j . 
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If Ox 7^ 0, then z/g would participate in the oscillations of atmospheric neutrinos 
with amplitude sin^2^a.. In the Super K data the number of observed electron neu- 
trinos is consistent with = (see Fig. IT^ . whereas an enhancement would be 
expected if there were ^ Ve oscillations (due to the 2:1 ratio of to Ve in the 
flux). Furthermore, the zenith angle distribution of the SuperK muon sample is in- 
consistent with oscillations involving z/g. Figure ITHl shows the SuperK allowed region 
in sin^ 6x — 5rn?^ plane. There is also slightly greater than 2a evidence in the SuperK 
data of hadronic showers from r decays \VAO\ I132j . consistent with the hypothesis 
that the primary oscillation of atmospheric neutrinos is ^ Vr- 
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Figure 18: Allowed regions in sin^ 0^ — Sm^ plane for three-neutrino oscillations from 
SuperK atmospheric data. From Ref. |131j . 

Another consequence of 7^ is that reactor z/g fluxes should exhibit disap- 
pearance due to the leading oscillation when L/Ey > 40 m/MeV. Data from the 
CHOOZ reactor experiment |2H1 {L ~ 1000 m, ii^j, ~ 3 MeV) place upper limits on 
6x for the values of 6m'^ indicated by the atmospheric neutrino data; similar limits 
have been obtained from the Palo Verde reactor experiment ^] (see Fig. [T^ . For 
the 6ml values obtained from the SuperK collaboration's latest two-neutrino analy- 
sis (slightly lower than from previous analyses), bounds on are quite sensitive to 
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Sml- While for Sml = 2.0 x 10~^ eV^, the 95% C. L. upper bound on sin^ 26'^: is 
0.2, for = 1.3 x 10~^ eV^, the corresponding bound is 0.36. Thus, it has become 
necessary to specify the 6m'^ for which a bound on 9^ is quoted. 

The K2K experiment, in which z/^ with energies typically ~ 1.4 GeV are directed 
from KEK to SuperK {L = 250 km), has measured a survival probability con- 
sistent with the atmospheric neutrino results, P{i^^ u^) = 0.70^0^1 tl^- The 
K2K allowed region, from the number of events and the spectrum shape combined, is 
consistent with the allowed region from the atmospheric neutrino data (see Fig. I19|). 
An approximate two-neutrino analysis that adopted the SuperK allowed regions in 
Fig. Uni in combination with K2K data found Sml = i2.0toi) x lO'^ eV^ at the 
95% C. L. |133j . Imposing this 6ml range as a prior in a three-neutrino analysis of 
CHOOZ, KamLAND and pre-SNO salt phase solar data yielded sin^26'2; < 0.17 at 
the 95% C. L. after marginalizing over 6ml !133j . 
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Figure 19: Allowed regions in sin^ 9a — 6ml plane from K2K, compared with the 
allowed regions from the SuperK atmospheric data |13()j . 

The MINOS experiment expects to detect atmospheric neutrinos via the z/^ 
and z/^ charged-current reactions and to determine the signs of the resulting charged 
leptons with a magnetic field, thereby separately testing oscillations of — > and 
i'^. 
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Since upward-going atmospheric neutrinos traverse a large fraction of the Earth's 
diameter, matter effects could be relevant. The dominant oscillation of atmospheric 
neutrinos appears to be z/^ — > Ur] for two-neutrino z/^ Vr oscillations there would 
be no matter effects. However, probability conservation in three-neutrino oscillations 
involves large matter effects which change the — > Uj. oscillation probability to a 
small extent. 

Specific relationships between the changes in the oscillation phase and the mat- 
ter density can lead to an enhancement of the oscillation probability, analogous to 
the classical phenomenon of parametric resonance |134j . Also, constructive quan- 
tum mechanical interference between the probability amplitudes for different density 
layers can give total neutrino flavor conversion |135j . Such enhancement phenom- 
ena generally require passage of the neutrino through the Earth's core. Conditions 
for observing matter effects in the Earth's mantle and core in atmospheric neutrino 
experiments are discussed in Ref. 

6 Absolute neutrino mass 

Neutrino oscillations tell us nothing about the absolute scale of neutrino masses. The 
standard technique for probing the absolute mass is to study the end-point region of 
the electron spectrum in tritium beta-decay. The effect of a nonzero neutrino mass is 
to suppress and cut off the electron distribution at the highest energies. The effective 
neutrino mass that could be determined in beta-decay is |137j 

ml = Y,\ye^?^l. (44) 

The present limit from the Troitsk j3] and Mainz jB] experiments is < 2.2 eV at 2a. 
Future sensitivity down to = 0.35 eV is expected in the KATRIN experiment |138j . 
which will begin in 2007. 

The sum of neutrino masses S = ^ mjy, can be probed in cosmology by measuring 
ujy = S/(94 eV) using the large scale structure (LSS) of the universe [TO]. Relativistic 
neutrinos do not cluster on scales smaller than they can travel in a Hubble time. Thus, 
the power spectrum is suppressed on scales smaller than the horizon when the neutri- 
nos become nonrelativistic. (For eV neutrinos, this is the horizon at matter-radiation 
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equality). The effect is subtle. Lighter neutrinos freestream out of larger scales and 
cause the power spectrum suppression to begin at smaller wavenumbers |139j : 



assuming almost degenerate neutrinos. On the other hand, heavier neutrinos consti- 
tute a larger fraction of the matter budget and suppress power on smaller scales more 
strongly than lighter neutrinos 170]: 



The galaxy power spectrum is influenced by the sum of neutrino masses, even down 
to 0.1 eV [10]. Analysis of the 2dF Galaxy Redshift Survey (2dFGRS) data found 
a limit of S < 2.2 eV on the masses of degenerate neutrinos, or about 0.7 eV for 
each neutrino |14Uj. An improved limit of S < 0.71 eV was obtained by the WMAP 
collaboration in an analysis of CMB data in conjunction with the 2dFGRS and the 
Lyman alpha forest power spectrum data fTT]. Since there are questions about the 
treatment of the uncertainties |141j in the Lyman alpha forest data, the WMAP col- 
laboration performed an analysis without this data and find the bound is strengthened 
to S < 0.63 eV Jlj. With only CMB and 2dFGRS data, other analyses found the 
limit on the summed neutrino masses to be about 1 eV |142j . In connection with the 
above limits, it is interesting that an argument relying on anthropic selection con- 
cluded that S ~ 1 eV so that neutrinos cause a small but nonnegligible suppression 
of galaxy formation |143j . In the future, lensing measurements of galaxies and the 
CMB by large scale structure may also provide a sensitive probe of S |144j . 

The Z-burst mechanism |145j provides another astrophysical probe of the absolute 
scale of neutrino masses |146j . The Greisen-Zatsepin-Kuzmin (GZK) cutoff energy 
(~ 5 X 10^^ eV) jl47j expected in the cosmic ray spectrum is absent in data from the 
AGASA experiment jl48j . although the GZK cutoff may be respected by the HiRes 
data |149j . A possible explanation for ultra high energy cosmic ray events above 
the GZK cutoff is resonant annihilation of ultra high energy neutrinos on the cosmic 
neutrino background to produce Z bosons which decay in a burst of about 20 photons 
and 2 super-GZK nucleons. The average energy of the secondary nucleon arriving at 




(45) 




(46) 
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Earth after travelling a distance D (~ 50-100 Mpc) is |150j 

^ 1021 eV X 0.8°/6Mpc 

m./0.1eV • ^''^ 
The neutrino mass scale constrained by atmospheric data from below and by cos- 
mology from above is suitable to initiate ~ 10^° eV air-showers. Higher would 
lower thus precluding Z-bursts as an explanation of the super-GZK events. Lower 

would necessitate an unrealistically large neutrino flux at the resonant energy. 
The speculative assumption of the Z-burst mechanism is that a substantial cosmic 
neutrino flux exists at ~ 10^^ eV. If the existence of this flux is confirmed at teraton 
neutrino detectors |151j . Z-bursts must occur. 

All neutrino masses are linked to the lightest mass by the values of and 5m1 
determined by the neutrino oscillation studies (152j . If the scale of the lightest mass is 



small, then the heaviest mass is approximately ^J \ Sm^ \ ~ 0.05 eV and a neutrino mass 
hierarchy exists. Since the sign of is unknown, there are two possible hierarchies, 
as illustrated in Fig. EDI The mass hierarchy is an important discriminant of neutrino 
mass models. If the scale of the lightest mass is larger than 0.05 eV, then the neutrino 
masses are approximately degenerate. 
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Figure 20: The patterns of relative mass differences in normal (left) and inverted 
(right) neutrino mass hierarchies. 

In the Standard Model with massive neutrinos and no other new physics, neutri- 
noless double-beta decay (Oz//5/5) probes the absolute mass, provided that neutrinos 
are Majorana particles; see Fig. |^ Numerous theoretical analyses have been made 
of what can be learned about the neutrino sector from Qvj3j3 |153| 11541 1155| . The 
decay rate depends on the v^-Ve element of the mass matrix jHH] : 

Mee=|^K>» ■ (48) 
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Figure 21: Neutrinoless double-beta decay mediated by Majorana neutrinos. 

The prediction is insensitive to 6^ and 6m1 because they are small. Setting 6^ = = 
5m^, the following relation between Mee and S is obtained for both hierarchies |154j : 

M,e = (2S - v/S2 + 35m2) jc^ + s^e^<^| /3 , (49) 

where is a Majorana phase. For a given measured value of M^e both upper (since 
6s 7^ 7r/4) and lower bounds are implied for E. These bounds are displayed in Fig. 1221 
The present upper limit on Mee is 0.35 eV at the 90% C. L. |157j . with an overall 
factor of 3 uncertainty associated with the Oi>f](3 nuclear matrix elements !15,8J . A 
detection of neutrinoless double beta decay, corresponding to Mee = 0.39 eV, has 
been reported |159j . but this experimental result is highly controversial |160j . 
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Figure 22: S vs. Mee for the normal (shaded) and inverted (dotted) hierarchies. For 
the inverted hierarchy, Mee > (Here, \6ml\ was taken to be 3 x 10"'^ eV^). 

The 95% C. L. bounds from tritium (3 decay and cosmology are shown. Adapted from 
Ref. [THij . 
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An extensive review of past and proposed Ou/SP experiments has been made in 
Ref. UnH]. Future experiments include CUORE {^^'^Te) [M], EXO (}^^Xe) [HSj . 
XMASS (i36Xe) [nOl, GENIUS (^^Ge) [HSI, Majorana (™Ge) HHSI and MOON 
(^^''Mo) jll8j . The upcoming experiments are expected to have sensitivity better 
than 50 meV, which is the critical mass scale of a/|5?7i^|. 

There has been speculation about detecting CP violation using OujSlS |155j . How- 
ever, it has been shown in Ref. |156j that this is impossible in the foreseeable future 
since there is no reliable method for estimating the uncertainty in the nuclear matrix 
elements. Moreover, the further the solar amplitude is constrained away from unity, 
the more stringent will be the precision requirement on the matrix elements for such 
a detection to be made even in principle |156j . Even under extremely optimistic as- 
sumptions, at best it may be possible to determine whether is closer to or to vr, 
corresponding to CP conservation. 

7 Supernova neutrinos 

Stars with masses above eight solar masses undergo collapse. Once the core of the 
star becomes constituted primarily of iron, further compression of the core does not 
ignite nuclear fusion and the star is unable to thermodynamically support its outer 
envelope. As the surrounding matter falls inward under gravity, the temperature of 
the core rises and iron dissociates into a particles and nucleons. Electron capture on 
protons becomes heavily favored and electron neutrinos are produced as the core gets 
neutronized (a process known as neutronization) . When the core reaches densities 
above 10^^ g/cm^, neutrinos become trapped (in the so-called neutrinosphere) . The 
collapse continues until 3 — 4 times nuclear density is reached, after which the inner 
core rebounds, sending a shock-wave across the outer core and into the mantle; see 
Fig. ESI This shock-wave loses energy as it heats the matter it traverses and incites 
further electron capture on the free protons left in the wake of the shock. During the 
few milliseconds in which the shock-wave travels from the inner core to the neutri- 
nosphere, electron neutrinos are released in a pulse. This neutronization burst carries 
away approximately 10^^ ergs of energy. However, 99% of the binding energy Eb, of 
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the protoneutron star is released in the following ~ 10 seconds primarily via /?-decay 
(providing a source of electron antineutrinos) , v^p^ |l(j5j and e^e~ annihilation and 
nucleon bremsstrahlung |166j (sources for all flavors of neutrinos including z/^, i/^, 
and z/^), in addition to electron capture. The neutrinos following the neutronization 
burst are the ones of interest in the following discussion. 




Figure 23: Schematic illustration of a supernova explosion. The dense Fe core col- 
lapses in a fraction of a second and gets neutronized (lower-left). The inner core 
rebounds and gives rise to a shock- wave (lower-right). The protoneutron star cools 
by the emission of neutrinos. From Ref. jl(i4j . 

We focus on charged current and Vf. interactions. We therefore cannot dis- 
tinguish between the different nonelectron species, and denote them collectively as 

[x = fi,T,p,,fY. The various cooling processes result in a state of approximate 
(within a factor of two or so) equipartition of energy with the luminosity of the elec- 
tron neutrinos L^^ being up to 10% larger than the that of the electron antineutrinos 

^However, note that weak magnetism can cause the luminosities and temperatures of the v^^t 
and Dfj_,r to differ by about 10% |1K7| . Practically speaking, this effect is too small to be detectable. 
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and 50 — 100% larger than that of |165j . The degree to which equipartition is 
violated can be quantified through constants and (3^^ which are defined by 

Lu, = P.M. = PuM. , (50) 
where 1 < ^1-1 and 1 < P^^: ^ 2. Perfect equipartition corresponds to = 
P.. = 1- 

Since the protoneutron star is opaque to neutrinos, it takes a few tens of sec- 
onds for them to diffuse out. The z/g and Ue interact with nuclear matter via both 
charged and neutral current reactions (with a smaller cross-section for Ue), while the 
z/j. experience only neutral current scattering. Consequently, the different species have 
neutrinospheres such that their radii obey R,y^ > Rp^ > R^^. Each neutrino species 
decouples at a temperature characterized by the temperature at the surface of its 
neutrinosphere. Based on this simple argument, which relies only on the well-known 
interaction strengths of neutrinos with matter, the hierarchy of average energies is 

(E.J < (EpJ < (E.J . (51) 

Early supernova (SN) models typically predicted jlti8j 
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The inclusion of additional energy transfer processes in modern SN codes indicates 
that the hierarchy of average energies is likely smaller than originally expected. There 
is evidence that nuclear recoils can lower (E.^) by as much as 20% |169j . Also, nucleon 
bremsstrahlung softened the spectra in the simulations of Ref. |170j . The cumulative 
effect is that {E^J / {EpJ is typically expected to be 1.1, and unlikely to be greater 
than 1.2 

The energy spectra of neutrinos can be modeled by pinched Fermi-Dirac distribu- 
tions. The unoscillated differential fiux (fiux per unit energy) at a distance D can be 
written as 

F = ^ f53) 

" 2A7TD^T^\LU{-e^'^)\ e^l^^-^c. + 1 ' ^^^^ 
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where a = i/g, z/g, z/^,, Lin{z) is the polylogarithm function and r/^ is the degeneracy 
parameter. The temperature of the neutrinos, T^, is related to {Ea) via {Ea) = 

Because of the extremely high density of the matter in the neutrino production 
region of a SN, all flavors of neutrinos start out in pure mass eigenstates. As the 
neutrinos stream out from the production region, they pass through a density profile 
that is well-represented by Vo(-R0/r)^ |172j . where Rq is the solar radius and Vq 
is a constant. Due to the wider range of densities that the neutrinos encounter, 
both the solar and atmospheric scales contribute to the oscillation dynamics. The 
hierarchical nature of the two scales (|5m^| <^ l^'^al) and the smallness of the mixing 
parameter sin^ 29^ imply that the dynamics can be approximately factored so that 
oscillations are governed by and sin^ 29^ at high densities (10'^ — 10^ g/cm^), and 
by 6m1 and sin^ 2^,, at low densities (~ 20 g/cm^ for the LMA solution) |173j : see 
Fig. 123 Transitions in the latter region are adiabatic. In the high density region, 
neutrinos (antineutrinos) pass through a resonance if > < 0). The 

jumping probability is the same for both neutrinos and antineutrinos |174j and is of 
the form Ph ~ ^-sin^ e^i\Srnl\/E^)^/^Vg^^ jl75j . Note the exponential dependence of Ph 
on sin^ 9^. That SN neutrinos provide a handle on the sign of (5m^ can be seen as 
follows. 



Figure 24: Schematic level-crossing diagram for neutrinos emitted by a SN in the case 
of a normal mass hierarchy, z/^/ and z/,-/ are basis states which diagonalize the (z/^, z/^) 
submatrix of the Hamiltonian governing the neutrino evolution. 
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For the normal and inverted hierarchies, the survival probability of electron an- 
tineutrinos is given by |173j 

P = Ae, (54) 

and 

p = PHPie + il-PH)sm^e,, (55) 

respectively. Here, Pie = Pie{E^, 6m'j., sin^ 26'^) is the probability that an antineutrino 
reaching the Earth in the ui mass eigenstate will interact with the detector as a Ue- 
If sin^ 29x <^ 10^'^, then P^ ~ 1 and the survival probabilities for the two hierarchies 
are the same. Thus, the normal and inverted mass hierarchies are indistinguishable 
for sin^ 29^ <^ 10~^. If sin^ 29^. > 10~^, for the inverted hierarchy p ^ sin^ 9^. < 0.05 
and the original electron antineutrinos have all been swapped for the more energetic 
fi and r antineutrinos by the time they exit the supernova envelope, resulting in a 
harder incident spectrum. Thus, the initial Ue spectrum would have to be softer for 
the inverted hierarchy than for the normal hierarchy. 

The detection of neutrinos from SN 1987A was momentous. The 11 events at 
Kamiokande II |176j and 8 events at the Irvine Michigan Brookhaven |177j detectors 
have lent strong support to the generic model of core collapse supernovae |178j . The 
significance of these few events provides a tantalizing glimpse into the physics poten- 
tial offered by a future galactic SN event. Despite the fact that only a few galactic SN 
are expected per century, the potential payoff is so huge that experiments dedicated 
to SN neutrino detection have been proposed [UB] . 

Attempts have been made to extract neutrino oscillation parameters from the 19 
SN 1987A events (see Refs. |179tll80[ll81j for recent analyses). However, conclusions 
drawn from these analyses depend crucially on the assumed neutrino temperatures 
and spectra. For example, it was claimed that the data favor the normal hierarchy 
over the inverted hierarchy provided sin^ 6'^; > 10~^ |18()j . but this conclusion was 
contradicted |181j . Table |21 shows the results of two-parameter fits in Ef,, and Tp^, for 
the normal hierarchy and the inverted hierarchy (with sin^ 26*2; = 0.01). In all cases 
the likelihoods are comparable, so the data do not favor one neutrino mass hierarchy 
over the other. While SN 1987A was of great astrophysical significance, it did not 
allow firm deductions about neutrino mixing. 
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Table 2: Best fit values for Eh and Tp^ from two-parameter fits to all the KII and 
1MB data. Results are for the cases of no oscillations, the inverted hierarchy with 
sin^2^2^ = 0.01 and the normal hierarchy. From Ref. jlSlj . 





Ef, (10^3 g^gs) 


Tp^ (MeV) 


If^i^^rnax ) 


no oscillations 


3.2 


3.6 


-42.0 


{E,J/{EpJ = 1.25 , sin2 2^, = 0.01 


3.1 


2.9 


-42.0 


{E^J/{Ei,^) = 1.25 , normal 


3.2 


3.4 


-41.9 


{E,J/{Ei,J = 1.4 , sin2 2^, = 0.01 


3.1 


2.6 


-42.0 


{E^J/{EpJ = 1.4 , normal 


3.4 


3.2 


-41.6 


{E,J/{EpJ = 1.7 , sin2 2^, = 0.01 


3.2 


2.1 


-42.0 


{E^J/{EpJ = 1.7 , normal 


4.2 


2.7 


-41.2 


{E,J/{EpJ = 2.0 , sin2 2^, = 0.01 


3.2 


1.8 


-42.0 


{E,J/{EpJ = 2.0 , normal 


5.8 


2.2 


-40.6 



Neutrinos from a galactic SN incident on a large water or heavy water detector 
could in principle provide much information on neutrino oscillations. A determination 
of 6x and the neutrino mass hierarchy from SN neutrinos is special in that degeneracies 
arising from the unknown CP phase 6 and whether 6a is above or below it /A do not 
contaminate it, i.e., the eight-fold parameter degeneracies that are inherent in long 
baseline experiments [UT] are absent. This cleanness results because (i) nonelectron 
fiuxes do not depend on the CP phase 6 |182j . and so SN neutrinos directly probe 
dx, and (ii) whether 6a is above or below 7r/4 is immaterial since this parameter does 
not affect the oscillation dynamics. 

Investigations of the effect of neutrino oscillations on SN neutrinos have been made 
in Refs. |173| 11831 1184j . Whether or not the mass hierarchy can be determined and 6^ 
be constrained depends strongly on how much {E^J / {EpJ is greater than unity |173j . 
The higher the value of (E^,^) / {EpJ, the better the possible determinations. As noted 
earlier, {E^J / {EpJ is expected to be about 1.1, and no larger than 1.2 [171j. Then, 
assuming that the rja values predicted by SN models are accurate, a safe deduction 
is that observations of a galactic SN at SuperK or Hyper-Kamiokande (HyperK) 
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could place either a lower or upper bound on Or^ if the neutrino mass hierarchy is 
inverted |184j . The hierarchy can be determined if it is inverted and sin^ 29^ > 10~^; 
for sin^2^2,' ^ 10~^, the survival probabilities of the electron antineutrinos are similar 
for both hierarchies rendering them indistinguishable even in principle. On the other 
hand, if the hierarchy is normal, neither can 6^ be constrained nor can the hierarchy 
be determined [184'. Nonetheless, the importance of the detection of a galactic SN 
should not be understated because of its major impact on the understanding of the 
SN explosion mechanism. 

Broadly speaking, the binding energy of the star and the temperatures of the 
different flavors of neutrinos and antineutrinos are determinable; see Fig. |2H1 |184j . 
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Figure 25: Determination of the binding energy E},, the supernova neutrino mean 
energies (temperatures) and sin^ 29^ for the normal mass hierarchy. The left-hand 
and right-hand panels correspond to data simulated at sin^ 29^ = 10~^ and sin^ 29^ = 
10~^, respectively. The cross-hairs mark the theoretical inputs, and the 90% and 99% 
C. L. regions are light and dark shadings, respectively. {E,y^) is the mean energy of 
the nonelectron neutrinos. From Ref. |184j . 
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Precise measurements of the temperatures can provide an unique window into the 
dominant microphysics of neutrino transport in addition to confirming the expectation 
that Tp^ < Ty^. With very large detectors it may even be possible to estimate how 
equipartitioned the energy is among the neutrino flavors. This knowledge would help 
reflne SN codes that predict different degrees to which equipartitioning is violated. For 
example, in Ref. |185j an almost perfect equipartitioning is obtained while according 
to Refs. |165t I186j . equipartitioning holds only to within a factor of 2. 

Since it is difficult to determine 9^ and the mass hierarchy simultaneously from 
galactic SN data, an interesting prospect is to consider what can be learned from a 
future SN if sin^ 2^^ is already known to be larger than 0.01 from a future reactor 
or accelerator neutrino experiment. With this information, SuperK would easily 
discriminate between the two hierarchies if the hierarchy is inverted. If the hierarchy 
is normal. Earth-matter effects on the SN neutrino flux could be used to verify that 
this is the case, with either a high-statistics detector like the proposed half-megaton 
HyperK or at a high- resolution scintillation detector jl87j . 

8 Model building 

8.1 Patterns of neutrino masses and mixings 

One of the important challenges in particle physics is to understand the spectrum of 
fermion masses. The mixing matrix in the quark sector, Vckm-, is given by the product 
V^Vdi where Vu and Vd are the unitary transformations applied to the left-handed up 
and down quarks to diagonalize the up and down quark mass matrices. Similarly, 
the mixing matrix that enters into neutrino oscillations is Vmns = ^lK/, where Vl 
and Vu are the unitary transformations applied to the left-handed charged leptons 
and neutrinos to diagonalize the charged lepton and neutrino mass matrices. In the 
quark sector, all mixing angles in Vckm are small and there is a mass hierarchy 
among the generations, whereas in the lepton sector (although not necessarily in 
the neutrino sector) a mass hierarchy exists with two large mixing angles and one 
small mixing angle in Vmns- A remarkable property of neutrino masses is that they 
are so much lighter than the charged leptons. Any theory of fermion mass must 
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reconcile the extreme differences between quark and lepton masses and mixings. A 
plethora of papers has addressed the problem of neutrino mass over the years in the 
context of different models, which we are unable to exhaustively cover here; for recent 
comprehensive reviews and more references, see Refs. |7H1 17^ 175] . 

Since absolute neutrino masses are not yet known, there are three possible mass 
patterns for neutrinos: (i) normal hierarchy (mi ^ 1712 <^ rn^), (ii) inverted hierarchy 
(m2 ^ mi ms), and (iii) quasi-degenerate (mi ~ m2 — ms). Because Vmns is a 
product of the mixing matrices for the charged leptons and neutrinos, the observed 
mixing in neutrino oscillations can originate from Vl, V^, or a combination of the 
two. Viable models exist with different combinations of mass pattern and origins of 
the mixing angles. 

In models where where the charged lepton mixing matrix is approximately diago- 
nal, Vmns derives directly from V^. If there are three Majorana neutrinos, then there 
are nine independent parameters in the mass matrix: three absolute masses and six 
mixing matrix parameters (see Eq. E)). Six of these may be measured in neutrino 
oscillations (three mixing angles, the Dirac phase, and two mass-squared differences). 
The absolute mass scale may be determined by measuring tritium beta-decay, Oz//?/5 
decay or the suppression of the matter power spectrum. The magnitude of the u^. — v^. 
mass matrix element (which depends on the three mixing angles, absolute masses, 
and two Majorana phases) may be determined from Oi//3/3 experiments (although the 
value of the associated nuclear matrix elements makes this measurement less than 
precise). Therefore the complete 3x3 mass matrix for Majorana neutrinos can- 
not be fully determined by experiment in the near future, and in practice neither 
of the Majorana phases will be well-measured. In models where the charged lepton 
mass matrix is nondiagonal, there are even more independent parameters. However, 
symmetry arguments can help to reduce the number of parameters. 

Soon after the initial SuperK discovery that atmospheric neutrinos oscillate with 
maximal or nearly maximal amplitude, it was noted that the neutrino sector might 
exhibit bimaximal mixing |188j . z.e., maximal or nearly maximal mixing of solar and 
atmospheric neutrinos. Then the mixing matrix has the unique form (up to state 
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redefinitions) 



V2 V2 

_ 1 1 

2 2 



\ 



Vmns 



1 



(56) 




2 V2 / 



Perturbations on this basic form can yield mixing tliat is not quite maximal, and 
can make Ves (= Sxe~'^^) nonzero. More recently, solar neutrino and KamLAND data 
disfavor maximal mixing in the solar sector (sin^ 26s < 0.95 at the 3cr level), and the 
emphasis is now on finding models that can give bi-large mixing. 

8.2 The seesaw mechanism 

A popular model for understanding the smallness of neutrino masses is the see- 
saw mechanism |189j . which is particularly well- motivated in a grand unified theory 
(GUT). In the one-generation version, the neutrino mass matrix in the vl-^r basis is 



where m^) is a Dirac mass and a right-handed Majorana mass. The eigenmasses 
are then approximately — m|,/m/j and tjir (the negative value of the lighter state 
can be made positive by a redefinition of the phases of the neutrino fields). If the 



the r lepton mass then niR ~ 10 GeV. Other interesting possibilities for mr, are 
the electroweak vacuum expectation value or the top quark mass which would imply 
rriR = 10^^ GeV (close to the GUT scale). Since heavy right-handed neutrinos exist 
in most grand unified models, a GUT/seesaw model is very attractive. 

In practice, there are three generations of neutrinos and the neutrino mass matrix 
for the light neutrinos is —MJ^M^^Md, where Mjj is the 3x3 matrix describing 
the Dirac neutrino masses (presumably related by symmetries to the charged lepton 
or quark masses) and Mr is the 3x3 matrix describing the right-handed Majorana 
neutrino masses. If Md has a hierarchical form (similar to the charged lepton mass 
spectrum), then the neutrino mixing angles in Vmns tend to be small unless there 
is an unnatural conspiracy between and Mr |19m 1191] . The choice of particular 
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heaviest light neutrino has mass of order ~ 0.05 eV, and the Dirac mass is 
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forms for Ml, and M^i can avoid this problem |19H 1192] . A more detailed discussion 
of such models and relevant references can be found in Ref. ^73] . 

8.3 GUT models 

Since GUT models relate quarks and leptons, the fact that quark mixing angles are 
large and lepton mixing angles are small is potentially a problem for these models. 
In fact, many GUT models predicted small solar neutrino mixing, now ruled out by 
data from SNO and KamLAND. However, there is a way around this difficulty, due 
to the fact that GUT theories relate leptons and quarks of opposite chiralities, e.g., 
right-handed down quarks are in the same fermion multiplets as left-handed charged 
leptons, and vice versa. Therefore the mixing we observe in the lepton sector is 
connected to the right-handed quark mixing, which is unknown and not constrained. 
So-called lopsided models |193j take advantage of this fact. The three-generation 
Dirac mass matrices for the charged leptons and down quarks have the forms 
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where e ^ cr ~ 1. The entries involving only the second and third generations 
come from specific Higgs Yukawa interactions, with no contribution to the middle 
diagonal. The entries with an "x" involve the first generation and are very small, 
and in many (although not all) models are generated by Froggatt-Nielsen diagrams 
mediated by exotic vector-like matter fields |194j : see Fig. 1211 The up-quark and 
neutrino mass matrices are approximately diagonal. Then the quark mixing element 
is small, but the leptonic mixing element V^3, which relates to the mixing of 
atmospheric neutrinos, is large. In lopsided models, the large atmospheric neutrino 
mixing comes from the diagonalization of the charged lepton mass matrix; the solar 
neutrino mixing angle arises from the structure of the three-generation right-handed 
Majorana neutrino mass matrix, also determined in some models by Froggatt-Nielsen 
diagrams. Most lopsided models are embedded in a GUT |195j . but some are not |19tij . 
In many cases, horizontal (z.e., family) sjTumetries determine the textures of the mass 
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matrices. A Monte Carlo study suggests that lopsided textures are favored by the 
data |197j . In any GUT framework, proper comparison with data can only be made 
after allowing for the renormalization group running of the neutrino mass terms in 
the Lagrangian |198j : for recent discussions, see Ref. |199j . Of particular importance 
are how zeroes in the mass matrices behave under renormalization and the stability of 
the large mixing angles |2UUj . For a more complete discussion and further references 
on GUT models, see Refs. [THll^. 
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Figure 26: Froggatt-Nielsen diagrams. Here a and h are family-indices and (x, X) ^''^^ 
vector-like fields of mass M and {9) is the vacuum expectation value of the flavor 
Higgses or flavons. The tree-level diagram (a) generates the mass of the third family 
and the lighter masses are obtained by 0{{6) /M)"- suppressions from diagrams (b) 
and (c). From Ref. [7U] . 

Lopsided models can yield either large or small solar neutrino mixing |i202] : for 
a discussion of which models can naturally yield the LMA solution, see Ref. |2U3j . 
The known SO(IO) models that satisfy all experimental data favor a normal hierar- 
chy |2Ulj . i.e., rris > mi, m2 in the neutrino sector, and therefore (5m^ > 0. Predictions 
for Ve3 vary [73j. 

Another interesting possibility is to assume that the unification group is repli- 
cated at the Planck scale, i.e., there is one copy for each generation. This leads to 
family- dependent U{1) symmetries as the theory breaks down to become the Standard 
Model at low energies. Some consequences of such models are discussed in Ref. |2M | . 
Models that utilize family symmetries can have either (i) both large mixing angles 
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in the neutrino sector deriving from the neutrino mass matrix |205j . or (ii) the large 
atmospheric neutrino mixing angle deriving from the charged lepton mass matrix and 
the large solar neutrino mixing angle deriving from the neutrino mass matrix ;206| . 
SUSY GUT models with i?-parity violation are discussed in Ref. |207j . 



8.4 Non-GUT models 

There are many alternatives to explicit GUT models (although it is often assumed 
that they could emerge from an unspecified grand unified theory). Some popular 
possibilities include: (i) the Zee model, (ii) models with low-energy new physics 
in which neutrino masses are generated as loops, such as supersymmetry with R- 
parity violation, (iii) so-called "democratic" models, (iv) models with triplet Higgs 
bosons, (v) models with specific textures for the neutrino mass matrix or special 
relationships involving the entries in the mass matrix, (vi) models with dynamical 
electroweak symmetry breaking, and (vii) models with large extra dimensions. In the 
non-GUT models, Mj, is unrelated to the charged fermion mass matrices, and hence 
a priori there are few constraints on its structure. In many cases, horizontal (family) 
symmetries can be used to provide constraints and produce a phenomenologically 
compelling model. 

Zee model 

In the Zee model |208j . which invokes radiative neutrino masses via a charged SU {2)i 
singlet and a Higgs field in the loop (see Fig. I77j) . the neutrino mass matrix has the 
approximate form 

/ A B \ 



(59) 



A 
\B J 

where A B. In Zee- type models, the diagonal elements of My are zero, and the 
remaining off-diagonal elements may be nonzero (but small compared to A and B). 
Such a texture can also result from an approximate L^. ~ L^^ — Lj. symmetry |2U9j : 
for some other possibilities see Ref. |21Uj . The mass matrix in Eq. (jS^ yields large 
mixing for both solar and atmospheric neutrinos (although some specific models yield 
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vacuum solar neutrino oscillations that are now excluded). The mass hierarchy for the 
Zee-type mass matrix is inverted, i.e., <^ mi,m2 and (5m^ < 0. The prediction 
of the Zee model that the solar neutrino mixing is nearly maximal mixing is now 
excluded by data. 



Figure 27: Neutrino mass is generated at one loop in the Zee model. 
New physics at low energy 

In models with low-energy new physics, neutrinos couple to a heavy fermion in the 
theory. Mass terms for the light neutrinos are generated by loop diagrams involv- 
ing the neutrino and the heavy fermion. If the heavy fermion coupling to the sec- 
ond and third generation neutrinos is larger than to the first generation, a normal 
mass hierarchy and large mixing for atmospheric neutrinos results. In the minimal 
supersymmetric extension of the Standard Model (MSSM) radiative neutrino mass 
generation is a direct consequence of R-parity violation |211j (see Fig. EH))- For spe- 
cific realizations to explain the neutrino anomalies see Ref. |212j . R-parity violating 
SUSY models that reproduce the neutrino mass and mixing parameters can have 
specific signatures in future collider experiments, such as lepton-number violating 
final states f2l'6\ . neutralino decay within the detector |214t 1215] . neutralino decay 
branching ratios |215t I21()] . multi-6-jet events with an isolated charged lepton |217j . 



Figure 28: The dominant one- loop diagram which generates Majorana neutrino 
masses for left-handed neutrinos in i?-parity violating models. The coupling A' vio- 
lates lepton number as well as i?-parity. 
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and multi-lepton events |218j . It is also possible that neutrino masses are generated 
at the two-loop level |219j . 



Flavor democracy 

In models with flavor democracy [220^ , the quark and charged lepton mass matrices 
have the form 

/l 1 l\ 



M oc 



1 1 1 
VI 1 1/ 



(60) 



while the neutrino mass matrix is approximately diagonal. This scenario also leads 
to large mixing for solar and atmospheric neutrinos and a normal mass hierarchy. 
The democratic model (and many other non-seesaw models) predicts that the solar 
neutrino mixing angle is maximal and the atmospheric neutrino mixing angle is large 
but not necessarily maximal, whereas the data indicate the opposite. 

Triplet Higgs bosons 

In models with triplet Higgs bosons, horizontal symmetries are used to constrain 
the texture of the neutrino mass matrix |221j . An example that uses an 5*2 x 5*2 
permutation symmetry in a four-neutrino theory (with one neutrino becoming heavy) 
is given in Ref. 



Textures or special relationships 

Since it unlikely that all nine parameters of the neutrino mass matrix can be de- 
termined, many studies have examined simpler structures with fewer independent 
parameters (with the assumption that the charged lepton mass matrix is diagonal). 
Some recent examples: 

(i) It has been shown that 3x3 Majorana mass matrices with three or more inde- 
pendent zero entries are excluded by current neutrino data but there are seven 
distinct textures with exactly two independent zeroes that are acceptable |223j . 
Note that since a Majorana mass matrix is symmetric, a reflected off-diagonal 



49 



zero is not counted as independent. Two of these textures lead to a normal 
mass hierarchy and the other five to a quasi-degenerate mass spectrum. In 
fact, it is possible to fully determine the neutrino mass spectra corresponding 
to these textures |224j . Several aspects of these seven matrices have been stud- 
ied in Ref. |225j . Some of these textures can be realized in a seesaw model 
with or without extra U{1) fiavor symmetries |22tij . These textures can also be 
obtained in models with three Higgs triplets and a sufficiently massive triplet 
Majoron fFTT. 

(ii) The weak-basis independent condition det(Mj, = 0) (which would be approx- 
imately true if the lightest neutrino was nearly massless) can also lead to a 
complete determination of the neutrino mass matrix |227j . 

(iii) Another possible condition that the neutrino mass matrix might obey is form 
invariance, UMU'^ = M, where f/ is a specific unitary matrix such that 
represents a well-defined discrete symmetry in the neutrino fiavor basis |228j . 
This condition leads to a variety of possible mass matrices, including all three 
of the allowed mass patterns j229j . If the discrete symmetry is the non-Abelian 
group y44, the mass pattern is quasi-degenerate |23Uj . 

(iv) If the sum of the neutrino masses is zero |231j . which can occur in models whose 
neutrino mass matrix can be expressed as the commutator of two matrices, 
only the inverted hierarchy and quasi-degenerate mass patterns are allowed by 
current neutrino data |232j . 

If the charged lepton mass matrix is not assumed to be diagonal, then there 
are more independent parameters. A simplifying ansatz may be used to reduce the 
number of parameters. For example, the Fritzsch ansatz |2Hl^j assumes Mn = M22 = 
Mi3 = M31 = for both the charged leptons and neutrinos, which can lead to 
acceptable phenomenology. The large mixing of atmospheric neutrinos can come 
from Vl EMI- 
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Dynamical electroweak symmetry breaking 

The seesaw mechanism can be reahzed in models with dynamical electroweak symme- 
try breaking (extended technicolor, or ETC, models). By suppressing the Dirac mass 
terms m^j in Eq. (j37j) . the heavy Majorana scale need not be so high, (typically the 
m£) terms can be much smaller than the ETC scale) |235j . Dynamical electroweak 
symmetry breaking due to a neutrino condensate has also been considered |236j . 

Extra dimensions 

Theories with large extra dimensions |237j have been postulated to avoid the hierarchy 
problem. In such theories, there is no very large scale {e.g., the GUT or Planck scale), 
and so the smallness of the neutrino masses cannot be obtained from a conventional 
seesaw mechanism; instead, it is a consequence of the suppressed coupling between the 
active neutrinos on the brane (the usual four- dimensional world) and sterile neutrinos 
in the bulk (Kaluza-Klein modes) or on other branes, associated with the small overlap 
of their wavef unctions. According to the particular model and coupling mechanisms, 
the neutrino masses can be either Dirac or Majorana. Some examples of models of 
neutrino mass in extra dimensions are given in Refs. |238j . However, no evidence of 
Kaluza-Klein modes, whose effects are like those of sterile neutrinos, has been found 
in the oscillation data. For a more complete discussion of theories with large extra 
dimensions, see Ref. [7^ . 

Alternatively, extra dimensions can be generated dynamically at low energies from 
a theory which is four- dimensional and renormalizable at high energies (a process 
called dimensional deconstruction |239j|). Acceptable neutrino phenomenology ap- 
pears to be possible in such a scenario |24Uj . 

Neutrino anarchy 

Finally, even though it is aesthetically pleasing to think that symmetries in one form 
or another account for the structure seen in the neutrino masses and mixings, it could 
be that an essentially random three-neutrino mass matrix can give the appropriate 
phenomenology |241j . In models with neutrino anarchy, it was originally thought 
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that large mixing angles are quite natural, and the value of 6^ could lie just below 
the current experimental bound. However, a recent study suggests that large mixing 
angles are not preferred if the mass matrix elements are truly random in a basis- 
independent way |242j . Statistical analyses of nonrandom structures have also been 
performed [TTITll^ . 

9 Leptogenesis 

The range of neutrino masses that the data suggest lends credibility to the seesaw 
mechanism. A direct consequence of the seesaw mechanism is leptogenesis |244j . 
This process generates a net lepton asymmetry = {rii — ni)/s because all of 
Sakharov's conditions |245j are met: (i) the heavy right-handed neutrinos Ni decay 
into a lepton-Higgs pair {IH) and into the CP conjugate pair with different partial 
widths, thereby violating lepton number; (ii) CP violation results from phases in the 
Yukawa couplings and neutrino mass matrices; (iii) the cosmological expansion yields 
the departure from thermal equilibrium. 

As the universe cools and the Ni drop out of equilibrium, their decays lead to a 
CP asymmetry |246j . 

r(iv, ^ iH) - T{N, ^ m*) 

T{Ni IH) + T{Ni ^ IH*) ■ ^ ' 

The lepton asymmetry generated will be Yi ~ J2^i/9*y where is the number of 
degrees of freedom. 

Since sphaleron j247j interactions preserve B — L but violate B and L |248j . the 
lepton asymmetry is partially converted to a baryon asymmetry Yb- In terms of the 
initial B -L jM?] . 

YB = aYB-L = ^YL, (62) 
a — 1 

where a depends on the processes in equilibrium. In the seesaw extended SM (MSSM), 
a = 28/79 (a = 8/23). Note that Eq. (jU^ is valid only for temperatures far above 
the weak scale (for a review see Ref. |25Uj ). 

An interesting aspect of leptogenesis is that it requires the low energy neutrino 
masses to be sufficiently light that the baryon asymmetry is not washed out by 
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neutrino- mediated L- violating scatterings |251j ^. This has been further explored in 
Refs. |253| I254j . Reference |254j finds that for the mass-squared differences relevant 
to solar and atmospheric neutrino oscillations, leptogenesis is the unique source of 
the baryon asymmetry provided the lightest of the heavy neutrinos is O{10^) times 
lighter than the other heavy neutrinos. 

Although a very appealing idea, leptogenesis is difficult to test [255]. The can 
be expressed independently of the neutrino mixing matrix V of Eq. Q |256j . Any 
connection between the CP phase 6 and the CP violation required for leptogenesis 
requires assumption about the texture of the Yukawa matrix, and is therefore model- 
dependent |257j . 

The only way to test leptogenesis directly is to constrain the Yukawa matrix and 
the masses of right-handed neutrinos. This can be achieved by searching for lepton 
flavor violating decays and the electric dipole moments of the charged leptons to 
which orders of magnitude improved sensitivity are expected in the near future |258j : 
fi — > 67, r /i7, r ^ 67 and /i 3e. Even so, any determination of the Yukawa 
matrix will depend on how precisely the Higgs sector is known and on assumptions 
about the GUT model. 

10 Future long-baseline experiments 

A summary of present knowledge of neutrino parameters is given in Table El along 

with the near future projects that will improve this knowledge. The mixing angle 6x, 

the Dirac CP phase 6, and the sign of Sm^ are as yet undetermined; their measurement 

will be the main goal of future long-baseline neutrino experiments. In this section, we 

discuss the the next-generation long-baseline experiments that are being considered 

after MINOS, ICARUS, and OPERA. 

^If leptons couple to a heavy SU(2)l triplet in addition to Ni, neutrino masses do not neces- 
sarily induce asymmetry washout effects and an upper bound on the neutrino masses cannot be 
placed 
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Table 3: Present knowledge of neutrino parameters and future ways of improving this 
knowledge. 



3-neutrino 

1 11 

observables 


Present knowledge 
(~ 95% C. L.) 


Near future 


Oa 


45° ± 10° 


P{v^. J^t^) MINOS, CNGS 


0s 


32.5° ±3.6° 


SNO NC, KamLAND 


dx 


<13° (for \bml\ = 2.0 x lO'^ eV^) 


P(z/e — > i^e) Reactor, P(z/^ —>■ Ve) LBL 


\5ml\ 


(2.0^^:^) X 10-3 eV2 


P{v^ v^) MINOS, CNGS 


sgn(5m^) 


unknown 


P{u^^u^),P{u^^u^) LBL 


\5m1\ 


{7.ltli) X 10-5 eV2 


P{ve -* v^) KamLAND 


sgn(5m^) 


+ (MSW) 


done 


5 


unknown 


P{u^^u,),P{u^^u,) LBL 


Major ana 


unknown 


Oz//?/? 


02 


unknown 


OZ//5/3 (if ~ 0,7r) 


03 


unknown 


hopeless 


nil, 


rriy < 1 eV 


LSS, OuPP, /3-decay 



10.1 Conventional neutrino beams and superbeams 

The near term agenda is to confirm atmospheric neutrino oscillations in accelerator 
experiments and improve the accuracy with which |5m^| and sin^ 29a are determined. 
Experiments that measure disappearance will establish the first oscillation mini- 
mum in P(z/^ — > i^fj.). The K2K experiment from KEK to SuperK ^H], a distance of 
L = 250 km, has begun taking data again following the restoration of the SuperK 
detector. The MINOS experiment from Fermilab to the Soudan mine |T7j, at a dis- 
tance of L = 730 km, will begin in 2005. It is expected to obtain 10% precision 
on \Sm'^\ and sin^ 26'^ in 3 years running. The CERN to Gran Sasso (CNGS) experi- 
ments, ICARUS [18j and OPERA also at a distance L = 730 km but with higher 
neutrino energy, are expected to begin in 2007. The appearance of Ur should be ob- 
served in the CNGS experiments, which would confirm that the primary oscillation 
of atmospheric neutrinos is z/^ — > i^^. 

The three parameters that are not determined by solar and atmospheric neutrino 
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experiments are Ox-, sgn(5m^), which fixes the hierarchy of neutrino masses, and the 
CP- violating phase 5. The appearance of v^. in z/g oscillations is the most 

critical measurement, since the probability is proportional to sin^2^a, in the leading 
oscillation, for which there is currently only an upper bound (0.2 for 5m\ = 2.0 x 10^'^ 
eV^ at the 95% C. L., from the CHOOZ reactor experiment [SHI)- combining 
ICARUS/MINOS/OPERA data, it may be possible to establish whether sin2 26'^ > 
0.01 at 95% C. L. ^25"9) . 

The study of z/^ — z/g oscillations also allows one to test for CP violation in the 
lepton sector j^. Intrinsic CP violation in the Standard Model requires both 6 
or vr and 9^ 7^ 0. In vacuum, the CP asymmetry in the channel, to leading 

order in the (5m^'s, is 

P{y^ u,) - P{u^ ^ z7e) ^ / sin2g,sin2g, \ / sin2AA ^.^^ 

For large-angle solar and atmospheric neutrino mixing the first factor on the right- 
hand side of Eq. (j63|) is of order unity. The existence of CP violation requires that 
the contribution of the sub-leading scale, A^, is nonnegligible, so large L/E^, values 
are essential. In practice, the CP conserving and CP violating contributions may 
have similar size |26()j . depending on the values oi L/ and O^- Furthermore, Earth- 
matter effects can induce fake CP violation, which must be folded into any deduction 
of 5 (on the other hand, matter effects are essential in determining the sign of 5m^). 

The standard proposed method for measuring CP violation is to compare event 
rates in two charge conjugate oscillations channels, such as —* z/g and z/^ z/g. 
However, there are three two-fold parameter degeneracies that are present when two 
such measurements are made, which may result in an overall eight-fold degeneracy (a 
parameter degeneracy occurs when two or more parameter sets are consistent with 
the same data): 

(i) The {6,6x) ambiguity W7\ [2611 1262( I263j . in which two different parameter 
pairs, (5, O^) and (5', 9'^), lead to the same values for P(z/^ z/g) and P(z/^ — > z/g). 

(ii) The sgn(5m^) ambiguity |^ 12611 1264t 1265] . where {S,9x) for one sign of 6ml 
gives the same values for the oscillation probabilities as (5', 6*^) with the opposite 
sign of 5m^. 



55 



(iii) The {6, 



I E 

a; 2 



9a) ambiguity I266j . where {6a, 5, Ox) gives the same values for 



the oscillation probabilities as (| — 9a, S', 9'^). This ambiguity exists because 



ambiguity vanishes at the experimentally preferred value for 9a (= f ). 

For each of these degeneracies, a duplicity in inferred values of 6 and 9x is possible; 
thus each of these degeneracies can confuse CP- violating parameter sets with CP- 
conserving ones, and vice versa. In many cases these degeneracies persist for all 
experimentally allowed values of 9x- An overview of these parameter degeneracies 
can be found in Refs. [HH 126 7j . 

There are two "magic" baselines that are valuable to resolve some of these param- 
eter degeneracies: 

(i) The detector is located at a distance that corresponds to the first peak of the 
leading oscillation (A^ = f ): 



Then the z/^ — Ue probability depends only on sin 5 and not cos 5 (see Eq. OHj) . 
the {6, 9,j.) degeneracy is broken, and 9^ is uniquely determined for a given 
sgn(5m^) and 9a |HH I261j . There is a residual (5, tt — 5) degeneracy, but this 
degeneracy does not mix CP violating and CP conserving solutions. Figure QUI 
shows the remaining degeneracies when LjEy is chosen to be at the first peak of 
the oscillation. Furthermore, if L is taken to be very long, large matter effects 
will break the sgn((5m^) ambiguity. The minimum distance needed depends on 
the size of 9^ and but generally L > 1000 km is required. 

(ii) The detector is located at a distance such that AAa = GpN^L/ \/2 ~ tt, which 
for the Earth's density profile implies L ~ 7600 km. Then only the leading 
oscillation term survives in Eqs. (j2Hl EHl ESI) , (IHUj). and the oscillation prob- 
abilities for Vf, appearance are independent of 5 and Sml jHI] . This allows an 
unambiguous measurement of 9x [268^ 1269] (modulo the {9a, | — 9a) ambiguity). 



the channel used to determine 9a, i^^ survival, only measures sin^ 29a- The 




(64) 
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For each of these magic basehnes, additional measurements at different L and/or 
values would be necessary to break the remaining degeneracies and determine the 
precise values of 5 and 6a- 



(a) Ev = 3.13 GeV, L= 1290 km 
^8 = 11/2 

sin^2e, = 0.01 



\5 = 0, n 



= 0.00298 



0.002 0.004 0.006 0.008 0.01 0.012 



6mJ = -3x10-''eV' 
sin^2e, = 0.0106 



T 

1^ 




5mg = 3x10 eV 
sin^2e^ = 0.01 

(b) Ev = 1 .77 GeV, L = 730 km 

0.002 0.004 0.006 0008 0.01 



(c) E 


v = 3.13 GeV, L = 


1290 km 




sin^2e^ 


sine^ 




0.01 

0.0055 


0.585 
0.811 




1/2 

*\ 1/2 








0,11 

^<x3it'2 

^3ll/2 



0.002 0.004 0.006 0.008 0.01 



Figure 29: Remaining degeneracies when Aq = | for (a) the (5, 6x) ambiguity, (b) the 
sgn(5m^) ambiguity, and (c) the {6a, f — 6a) ambiguity. In (a), each value of sin^26'a; 
describes a distinct line in probability space. In (b), the ambiguity in sin^ 26x is small, 
but in the overlap region there is still an ambiguity in sgn(5m^) and a corresponding 
large uncertainty in 5. In (c), the ambiguity in 5 is small, but there may be a large 
uncertainty in sin^ 26 x when 7^ f • In all cases there remains a (5, vr — 5) ambiguity 
since only sin 5 is being measured. Adapted from Ref. |61j . 
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Precision measurements of 5 and 6^ can be made at future off-axis jHll neutrino 
beam experiments proposed for the Main Injector at Fermilab (NuMI) 12701 \'271\ 
\'2T2\ and the Japan Hadron Facihty (JHF), also called the Japan Proton Accelerator 
Research Complex (J- PARC) jHOl- Compared to on-axis beams, off-axis beams have 
a much narrower energy spectrum, smaller beam contamination, and a suppression 
of the high-energy tail that results in lower backgrounds to u^. events in the detector. 
They also are well-suited for multiple detectors (a detector cluster |273j ) that allows 
more than one measurement to be made simultaneously. 

Strategies for the future include superbeams I274j . with upgrades of the neu- 
trino flux by a factor of 4 for SuperNuMI (SNuMI) |275j . a factor of 5 for Super JHF 
(SJHF) [nO], and a factor of 5 for a low-energy option for CNGS |276j . which will allow 
smaller values of 6^ to be probed. A high-intensity neutrino beam at Brookhaven is 
also being considered |277j . Many different superbeam scenarios have been consid- 
ered: 

(i) Having detectors at different off-axis angles allows one to modify both the base- 
line and neutrino energy. Multiple detectors utilizing an off-axis beam such 
as the one at Fermilab can shift degeneracies to sin^ 26'^^ < 0.01 — 0.001 |273j . 
Alternatively, having detectors at two different on-axis distances from the same 
superbeam can provide multiple measurements that help to remove parameter 
degeneracies. One such possibility is to combine a moderate distance such as 
JHF to SuperK {L = 295 km), with a much longer distance such as JHF to 
a detector near Beijing (L ~ 2100 km) |278j . One specific proposal is for a 
SJHF to HyperK experiment to run 2 years with neutrinos and 5 years with 
antineutrinos from, and then run with 5 years with neutrinos from SJHF to a 
very large Water Cherenkov detector near Beijing. This very ambitious scenario 
would require a separate beamline for the 2100 km measurement, but could re- 
solve the {6, Ox) ambiguity down to sin^ 26x = 0.005 |279j . Finally, one can vary 
the beam energy at the same baseline |273t 1280] . 

(ii) Two superbeam experiments do significantly better than one in parameter de- 
terminations |28H I282j . For example, the combination of data from a SJHF 
to SuperK experiment (6'ofr-axis = 2 deg, E^, ~ 0.6 GeV, L = 295 km, 22.5 kt 
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water Cherenkov) and a SNuMI to southwestern Ontario experiment (1 deg, 
El, = 1.8 GeV, L ~ 900 km, 20 kt low-Z calorimeter), both with 2 years z/^ 
running and 6 years z/^ running, would be sensitive to the sign of (5m^ and to 
CP violation for sin^ 2^^^ > 0.03 |281j . Running with only z/^ at SNuMI and 
SJHF may allow one to determine sgn((5m^) if 6^ is not too small |283j . although 
6x and 6 will not be well-measured without z/^ data. 

(iii) Another approach is to use a wide-band superbeam, for example from Brook- 
haven to a National Underground Science Laboratory |277j . The measurement 
of quasi-elastic events allows a determination of the neutrino energy with rea- 
sonable precision. The lower-energy events are more sensitive to the 6 terms 
in the oscillation probability, while the higher-energy events are more sensitive 
to the sign of Sm"^. Binning the quasi-elastic events is roughly equivalent to 
running many narrow-band beams simultaneously, which can help to resolve 
neutrino parameter degeneracies. In principle only a neutrino beam is required. 
However, running with an antineutrino beam would provide essential confirma- 
tion (especially of CP violation) and probe lower in 6^ (if Sml < 0). 

In any long-baseline neutrino experiment there is a trade-off between detector 
size and the ratio of the z/g CC signal events to background: generally speaking, 
detector technologies that allow a bigger reduction in background cannot be built 
as large [B^. Four types of detectors that have been studied for z/^ detection 
are (i) water Cherenkov (backgrounds of order 10~^ of the number of unoscillated 
CC events, maximum fiducial volume of order 500 kt) [Hniinnij (h) iron scintillator 
(3 X 10-3, 50 kt) j2nil2Hlll2Hni, (iii) liquid argon (3 x lO^^, 50 kt) ^ and (iv) low- 
Z calorimetric |272j . There is also an additional background of order 3 x 10^'^ due 
to z/g contamination in the beam. The larger water Cherenkov detectors generally 
do better when the neutrino flux is less (such as for a conventional beam before 
superbeam upgrade or for baselines > 4000 km where there is a large 1/ B? fall-off of 
the flux), whereas the smaller detectors that can measure positions on a finer scale 
generally do better when there is more flux (such as with superbeams or for baselines 
below 4000 km). 



59 



Another option being considered is an intense v^, beam (with neghgible contamina- 
tion) from radioactive ion decay leading to a so-called beta-beam |28(jj in a proposed 
experiment from CERN to Frejus with a baseline of 130 km. 

10.2 Future reactor experiments 

A different approach to determining without the complication of parameter degen- 
eracies is to measure P[pe ^e) at a reactor experiment using two larger versions 
of the CHOOZ detector, one near and one more distant (< 10 km) from the reac- 
tor [2HZI EHHl ESH! ■ Table m lists the detector distances for three proposals. Ignoring 
terms cubic or higher in the small parameters Qx and A^, the oscillation probability 
is given approximately by 

P(z/e ^ z/e) = 1 - sin^ 20^ sin^ - c\ sin^ IQ, sin^ A, . (65) 

It is independent of 5, Qa-, and sgn((5m^). Discovery down to sin^26'j; = 0.01 at 
90% C. L. may be possible |2871 12881 12891 for a reactor experiment with exposure 
of about 400 t-GW-yr, where the tonnage refers to the detector size and GW to the 
reactor power; an actual measurement of sin^ Id^ is possible for sin^ IQ^ — 0.05. A 
comparison of the sensitivity to sin^ 29x in superbeam and reactor experiments is 
shown in Fig. 1201 When a reactor measurement is combined with results from long- 
baseline experiments, it may also be possible to determine 6 and 9a. 

Table 4: Proposed reactor neutrino experiments for measuring 6^ using two detectors. 



Site 


Li (km) 


L2 (km) 


Krasnoyarsk |287j 


0.1 


1.0 


Kashiwazaki |288j 


0.3 


1.7 


Diablo Canyon j2H 


SI 0.15 


1.2 



10.3 Neutrino factories 

A neutrino factory (NuFact) [HH] is the ultimate technology for neutrino oscillation 
studies. Muons would be stored in a fiat oval ring. Their decays will give neutrino 
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Sensitivity to sin^29x 




10 ^ 10 ^ 10 ' 

sin"2fti sensitivity limit 



Figure 30: Sensitivity to sin^ 26x at 90% C.L. in reactor experiments with 400 t-GW- 
yr (Reactor-I) and 8000 t-GW-yr (Reactor-II), and for a JHF-SuperK experiment 
assuming 5ml = 7^ lO'^ eV^ (LMA-I) and 5ml = 1.4 x 10"^ eV^ (LMA-II). The left 
edges of the bars show the sensitivities assuming statistical uncertainties only, while 
the right edges of the bars show the sensitivities as the uncertainties from systematics, 
parameter correlations, and parameter degeneracies are progressively included. From 

Ref. pnnj . 

beams in the directions of the straight sections of the storage ring. Stored muons 
of energies 20 GeV and above are needed; energies as high as 50 GeV have been 
considered in design studies. A decaying /i+ in the ring yields both z/^ and Ue] detection 
of the charge of the final state lepton in a charge-current event allows one to determine 
the initial neutrino flavor. Table El shows the six oscillation channels possible with 
stored /i"*"; the six charge-conjugate channels can be tested using /i~ decays. 

The neutrino spectrum from a NuFact ranges from zero to the stored muon energy, 
Efj_, with a broad peak at O.T-E^ (O.G-E^) for muon (electron) neutrinos when the 
muons are not polarized. The neutrino flux in the forward direction is approximately 
no7^/(7rL^), where riQ is the number of decaying muons in the straight section of the 
ring, 7 = E^/m^, and L is the baseline. An entry-level NuFact produces a time 
integrated ~ 10^° and a high-performance NuFact has ~ 10^^. Early studies 
of the capabilities of a NuFact can be found in Refs. I291j . See also more recent 
study group reports in Ref. |292j and the reviews in Ref. |293j . 
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Table 5: Signals for oscillation channels assuming a decaying jj in a neutrino factory. 



Channel Detect Nomenclature 









right- 


sign fi survival 




Ve 


e~ 


right- 


-sign e appearance 




Vt 


T~ 


right- 


sign T appearance 




Ve 


e+ 


wronj 


5-sign e survival 


ye 


Vy. 


^+ 


wronj 


i;-sign /i appearance 






r+ 


wronj 


5-sign r appearance 



Golden channel: z/g ^'/^ 

Since the sign of the detected lepton is critical for determining the oscillation channel 
being observed, most studies have focused on the z/g and z/g — > oscillation 

channels with final state muon detection. Employing a magnetized iron detector 
allows one to determine the sign of the detected charged lepton, and the backgrounds 
are quite small, of order 3 x 10~^. This small background compared to Ve detection 
in a superbeam, plus the fact that NuFact neutrino fluxes can be one or two orders 
of magnitude greater than that of a superbeam, are the two main reasons a NuFact 
is clearly superior. 

There have been many many studies of the physics capabilities of a NuFact using 
muon appearance jHIl I265| I294| I295| I296| I297j (for a discussion of the physics that 
can be done with electron appearance, see Ref. j298l). By comparing Vf, and 
event rates, factoring in Earth-matter effects, very precise determinations 
of the oscillation parameters dx, S, and sgn(5m^) can be made. Figure 1^ shows the 
ratio of antineutrino to neutrino muon appearance events versus baseline for > 
and < for sin^ 29^ = 0.004 and several values of 6. The different signs of 
are clearly distinguishable when L > 2000 km, and these measurements are especially 
sensitive to the amount of CP violation when L ~ 3000 km. In practice, if a NuFact 
is run with fi^ about twice as long as with fi~^, then the total number of CC events in 
the neutrino and antineutrino channels will be about the same since cross section 
is about half that of the Ue cross section. 
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Figure 31: Ratio of antineutrino to neutrino appearance events versus baseline in a 
neutrino factory for sin^ 26^ = 0.004 and several values of 5. Both 5rn?^ > and 
< cases are shown. From Ref. |295j . 

One obstacle to making precise measurements of 5 and Ox and determining the 
sign of is that the other neutrino mass and mixing parameters, i.e., 6a, Os, ^'^s, 
and Sm^, may not be precisely known. Measurements of z/^ survival in a superbeam 
or NuFact will reduce current uncertainties in 6s and Sm^, while future KamLAND 
and SNO measurements will reduce the uncertainties in Og and especially 5m^. 

As with superbeams, there is the possibility of having an eight-fold parameter 
degeneracy using neutrino and antineutrino event rates in a NuFact. The sgn((5m^) 
ambiguity can be resolved by choosing a baseline > 2000 km. The ambiguity is easier 
to resolve for large 9^ (due to the larger matter effect) and small (5m^ (due to the 
smaller size of the CP violating term in the oscillation probability). There have 
been a number of different proposals for resolving the {6, Ox) ambiguity in the golden 
channel: 

(i) Since a NuFact has a broad spectrum of neutrino energies, measuring the energy 
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of the detected muon gives information about the modulation of the oscillation 
probability with E^. A 10% muon energy resolution is sufficient to remove 
the {S,6x) ambiguity |296j . A combined fit with the survival channel can 
also improve the measurement of 6, 6^, and sgn((5m^). An understanding of 
parameter correlations and degeneracies is essential for extracting meaningful 
constraints from the data |2ti31 1296] . 

(ii) Another idea is to combine Ue —>■ and —>■ i^ti measurements from neutrino 
factory experiments at two baselines. Having one detector at L ~ 3000 km and 
another at L ^ 7300 km would provide good discrimination between degener- 
ate solutions for a wide range of 6 and 9^ (see Fig. 1221) • Measurements at 
7300 km, near the magic baseline, where the 6 and 6m1 dependence is mini- 
mal jni], would provide an unambiguous measurement of the leading oscillation 
amplitude sin^ 9a sin^ 29^ |268j . With an oval design this scenario would require 
two separate runs for both stored fi~^ and fi~ , whereas a triangular configuration 
could allow data to be taken simultaneously at two baselines. 

(iii) Since a superbeam will most certainly be a precursor to a NuFact, it is quite 
natural to combine data from a superbeam and a NuFact to help resolve the 
{6, 9x) ambiguity. Studies show that it is possible to remove this ambiguity for 
s\.v?29x > 0.0005 |299j . In this sense, superbeams and neutrino factories are 
complement ar y. 

Some comparisons of superbeam versus NuFact performance are given in Refs. [HSl 
I2(j3j . It should be noted that even if = 0, subleading terms in the oscillation 
probability associated with 5m1 can lead to observable effects in appearance experi- 
ments |295j . 

A summary of typical capabilities of future long-baseline experiments is given in 
Tableini A neutrino factory would have sensitivity down to 5 x 10~^ and possibly lower 
in s\.v?29x for both sgn(5m^) and CP-violation determinations. A 1% determination 
of (5m^ should be possible at a NuFact. Thus, it appears that precision reconstruction 
of the neutrino mixing matrix would be possible with a neutrino factory even for very 
small values of 9^- 
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Figure 32: Fits to 6 and 0^{= Oi^) at a neutrino factory using hypothetical results 
from two baselines, 2810 km and 7332 km, for several input values of 6 and 6^ (and 
9a = vr/4). The three curves in each case represent the 90%, 95%, and 99% C. L. 
ranges of allowed parameters. Expected uncertainties in the oscillation parameters 
have been included, in addition to a 1% uncertainty in the matter density. From 
Ref. jnZl, in whose notation Aml^ > implies Sm^ > 0. 

Table 6: Approximate 3a reaches in sin^ 29x of future neutrino oscillation experiments. 

Reach in sin^ 29 x 



Experiment 


Discovery 


sgn(5m^) 


CP violation 


Reactor 


0.01 






Conventional u beam 


0.01 






superbeam 


0.003 


0.003 


0.03 


Entry-level NuFact 


0.0005 


0.0003 


0.002 


High-performance NuFact 


0.00005 


0.0001 


0.0005 



Silver channel: 

Another advantage of a NuFact over a superbeam is that it has available the oscillation 
channel Ve ^r- The oscillation probabilities for z/g and z/g can be 

obtained from those for z/g and z/g by the transformations sin 9a ^ cos 9a 
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and 6 —6. Comparing the muon and tau channels in a NuFact therefore allows one 
to (i) help resolve the {6a, ^ — 6a) ambiguity since the leading term in the oscillation 
probability is proportional to sin^ 6*^ for and cos^ 6*^ for — > [61 , and 

(ii) help resolve the (5, 6^) ambiguity since the 5 dependence is different in the two 
channels jSOOj . 

10.4 T and CPT symmetries 

If CPT is conserved, then CP violation implies T violation, which can be measured, 
e.g., by comparing z/g to z/g. There have been many phenomenological 

studies of T violation in the literature |267t I297| l3Ulj . Unlike CP violation, matter 
does not induce T violation in a long-baseline neutrino oscillation experiment, due to 
the symmetric matter distribution {i.e., the matter distribution is the same from the 
detector to the source as for the source to the detector), although matter can modify 
the amount of T violation. 

If CPT is not conserved, then P(z/q, — > up) is not necessarily equal to P{i^i3 — * i^a) 
in a vacuum; matter can induce fake CPT violation (for a study of matter-induced 
CPT violation see Ref. |3()2j ). In a NuFact, CPT violation can be tested down to 
very low levels by comparing the survival channels z/^ —>■ and z/^ — > u^. There 
have been a number of studies of possible tests of CPT violation in the neutrino 
sector |3U3j . 

11 The outlier: LSND 

The focus of this review so far has been on three-neutrino phenomenology. However, 
the results of the LSND experiment may cast some doubt on the three-neutrino 
picture. The LSND experiment found evidence for z/^ — ^ z/g oscillations at 3.3cr 
significance (oscillation probability (2.64±0.67±0.45) x 10~^ [3^] in data on fi~ decays 
at rest taken from 1993-1998. Evidence for z/^ — > z/g oscillations was found at lesser 
significance from 7r+ decay in flight, with oscillation probabilities (2.6±1.0±0.5) x 10~^ 
in the 1993-1995 data HOI and (1.0±1.6±0.4) x lO^^ ^j^g 1996-1998 data (see the last 
paper of Ref. ISH])- There was a significant difference in the analysis of the decay in 
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flight data in the two time periods due to changes in the neutrino production target, 
so the two Vf, samples were not combined. The KARMEN experiment [H] 

also searched for z/^ — > v^. oscillations with a null result; the KARMEN data rule 
out a large fraction of the LSND allowed region, but still allow a limited region of 
oscillation parameters l 3U4j . The Bugey reactor experiment |3U5j . which tests the 
oscillation channel v^. z/g, excludes the part of the LSND region with sin^ 20^ > 
0.04. In a two-neutrino parameter space, the indicated oscillation parameters from 
a combination of LSND and KARMEN data that are consistent with the constraint 
from Bugey are b'm\ ~ 0.2 - 1 eV^, sin^ 2^l ~ 0.003 - 0.04 and bm\ ~ 7 eV^, 
sin2 2^i ~ 0.004 at the 90% C. L. (see Fig. ESI)- The bulk of the allowed region is 
a narrow band in (sin^2^i, plane lying along the line described approximately 
by sin2 2^L(5mi)^-*^^ = 0.0025, between bm\ = 0.2 and 1 eV^. The MiniBooNE 
experiment jlH| will search for oscillations over the entire parameter space 

allowed by the LSND results. 




Figure 33: The 90% C. L. allowed region (shaded) from a combined fit to — z/g data 
from LSND i39| and KARMEN g^. The 90% C. L. allowed region from LSND alone 
is unshaded. The 90% C. L. exclusion regions from KARMEN, Bugey (z/g —>■ Vx) |305j . 
CCFR (z/^ ^ z/g) inOn] and NOMAD (z/^ ^ z/g) [Ml and the expected 90% C. L. 
sensitivity of the MiniBooNE experiment jlH] are also shown. From Ref. |304j . 
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11.1 Four neutrinos? 

The LSND parameters are very different from tlie oscillation parameters that explain 
the solar and atmospheric neutrino data; in particular, \Smj^\ ^ l^^al' I'^^^sl- Since 
a theory with three neutrinos has at most two independent mass-squared differences, 
a third 6171^ scale suggests that there may be a fourth light neutrino participating in 
neutrino oscillations Because neutrino counting experiments indicate that there 
are only three light active neutrinos a fourth light neutrino must be sterile, i.e., 
it does not participate in the weak interactions |3U8j . 

In Section 2, we saw that an extra fully thermalized neutrino is strongly disfavored 
by BBN |78|. The only way the LSND sterile neutrino can be reconciled with BBN 
is if its mixing with the active neutrinos occurs only after the active neutrinos have 
decoupled from the — 7 plasma. Then, the LSND neutrino is not thermalized and 
the effective number of neutrinos remains equal to 3. An electron neutrino asymmetry, 

Le = ^^^^^-^ ^ 0.01 - 0.1 , (66) 

accomplishes this |309j and simultaneously improves the agreement between the BBN 
prediction for the primordial ^He abundance and the observationally inferred value |310j 
Also note that since the production of the LSND neutrino is suppressed, cosmological 
bounds on ^ m,^ pertain only to the active neutrinos. A confirmation of the LSND 
signal by MiniBooNE could be interpreted as a hint for a large neutrino asymmetry 
in the universe. 

There are two types of mass spectra possible in four-neutrino models ^HIIS]- In 
3 + 1 models, one mass eigenstate is separated from a nearly degenerate triplet of 
mass eigenstates by 6m\] the triplet has a mass ordering like that of a three- neutrino 
model. The well-separated mass eigenstate can be either lighter or heavier than the 
other three, and the triplet can have a normal or inverted hierarchy, so there are four 
possible variations of 3 + 1 models. In 2 + 2 models there is one pair of closely-spaced 
mass eigenstates separated from another closely-spaced pair by 5m|; one pair has a 
mass-squared difference of Sm^ and the other Sm^, and the solar Sm^ can be in either 
the upper or lower pair. Figure OH shows the six possible mass spectra with four 
neutrinos. 
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Figure 34: The six possible mass spectra in four-neutrino models. 
11.2 Four neutrino models 

The 3 + 1 models are a straightforward extension of a three-neutrino model: the 
three active neutrinos have mass-squared differences and mixings similar to those in 
a three-neutrino model, and the sterile neutrino state has only small mixing with 
active neutrinos. However, 3 + 1 models have trouble accounting for the LSND 
results and simultaneously obeying the constraints of earlier accelerator and reactor 
experiments HI] . This can be demonstrated as follows. 

Assume a neutrino mass spectrum such that the nearly degenerate triplet of mass 
eigenstate is lighter than the remaining state and exhibits a normal hierarchy (the first 
spectrum shown in Fig. lMj) : similar conclusions can be drawn for the other three 3 + 1 

spectra. Then Sml^ ~ 6ml2 — Sml^ = 5m\ ^ <5^32 — ^"^ii = ^^'i ^ '^^21 — ^'^^si 
and the oscillation probabilities for the leading oscillation, due to 5m\, are 

A\V^A?\VeA? sin^ /\l , (67) 
1-4|V;4|'(1- |V;.4|')sin2Ai, (68) 
l-4|V;4r(l- |K4r)sin2Ai, (69) 

where A/, = 5m\L/{AEy), analogous to Eq. (0). At L/Ey values appropriate for 
atmospheric neutrinos, 

P{v^ ^ zy^) ^ 1 - 4|r^3|'(l - l^^sl' - KA'') sin^ A, , (70) 
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and for solar neutrinos 

P{iy, ^ z/e) ~ l-4|reinre2|'sin2A,. (71) 

There are very stringent limits on disappearance from the CCFR |306j . NO- 
MAD |3U7j and CDHS accelerator experiments that constrain |V^4| to be very 
small or very close to unity via Eq. ()68|) . However, if |V^4| is close to unity then 
the amplitude of atmospheric neutrino oscillations cannot be as large as required by 
observation (see Eq. EOl); therefore, |V^4p <^ 1. Similarly, the Bugey reactor ex- 
periment |3(J5j puts severe limits on disappearance that constrain |Ve4| to be very 
small or very close to unity, and the observation of large angle mixing in solar neutrino 
oscillations therefore implies \ Ve4\ cannot be close to unity (see Eg. ITTj) . so |Ve4p <^ 1. 

Since the oscillation amplitude for the LSND experiment is 4|V^4p|\44p, it has 
an upper limit of approximately one-fourth of the product of the CDHS and Bugey 
oscillation amplitude bounds (see Eqs. lUTHU^ . In practice, the limits on |V^4| and 
\Vei\ depend on (5m|, as does the allowed oscillation amplitude from LSND, so a 
comparison must be made for each value of (5m| . Early analyses |^ concluded 
that the 3 + 1 model could not consistently explain the LSND, accelerator, and 
reactor data. Upon the release of the final data analysis by the LSND collaboration 
(the last paper of Ref. [Sn], in which the central value of the average oscillation 
probability decreased from 0.0031 to 0.0026), it was thought that perhaps the 3 + 1 
models would be revived [IHllini- Subsequently, a Bayesian analysis of all relevant 
data was made |3T2] . Figure ESI shows the incompatibility of the combined accelerator 
and reactor upper limit on the oscillation amplitude in 3 + 1 models and the region 
allowed by LSND and KARMEN. An analysis of all relevant data yielded the best 
overall goodness of fit to be 5.6 x 10"'^ jl7j. Thus, the viability of 3 + 1 models is 
tenuous. 

The 2 + 2 models are not a simple extension of a three-neutrino model, since the 
removal of the sterile neutrino does not leave the standard three-neutrino mass spec- 
trum. If the lower pair of neutrino mass eigenstates are primarily responsible for solar 
neutrino oscillations, then we have 5m|2 — ^^32 — "^"^si — S'^li = Sm'j^ ^ 6m1^ = 
5m?^ ^ 5m\^ = dm^ (the case where the lower pair of mass eigenstates are primarily 
responsible for atmospheric neutrino oscillations leads to similar conclusions). Then 
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Figure 35: Upper bounds on the LSND amplitude in 3 + 1 models (4| V^4p| V"e4p) 
from accelerator and reactor data |312j . The shaded regions are allowed by LSND 
and KARMEN at the 95% C. L. jHHij . 

the oscillation probabilities for the leading oscillation are 

P{u^^ue) ^ 4\V,sV:, + V,,V:fsm'AL, (72) 
P{iy^^u^) ^ i_4(|r^3r + |V^M4r)(l-|^M3r-|V;4r)sin2Ai, (73) 
Piue^i^e) - l-4(|V;3p + |K4|')(l-|K3r-|V;4r)sin2Az., (74) 

At an L/ appropriate for atmospheric neutrinos, 

P{u^ ^ zy^) ~ 1 - 4|F^3n^M4l' sin^ A, , (75) 

and the oscillation probability for solar neutrinos in the 2 + 2 case is the same as in 
the 3 + 1 case fEa. r7T|) . 

If Ve is primarily connected to vi and V2i and is primarily connected to z/3 and z/4, 
I.e., \Vei\^ + |V;2p, \V^,-i\^ + |1/^4p ^ 1 > \Ve^?, \VeA\^ , \V^,i\^ , \V^2? , then it is possible 
to simultaneously fit the solar, atmospheric, and LSND data in 2 + 2 models [121 El- 
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For example, if |V^3p + |V^4p = 1 and + — 1, then it is easy to have 

large mixings of solar and atmospheric neutrinos while suppressing oscillations in 
the sensitivity regions of CDHS and Bugey. Since -P(t'^ — > i^e) and P{ue ^ ^e) are 
both proportional to the second power of small parameters (Ves and Vei)-, the LSND 
oscillation amplitude is approximately constrained by the Bugey bound in the 2 + 2 
case (in contrast to one-fourth of the product of the CDHS and Bugey bounds in the 
3 + 1 case). This less severe constraint is readily satisfied for a wide range of 5m\\ 
in fact, if either Ves or Ve4 is zero, then the constraint exactly reduces to the simple 
two-neutrino Bugey constraint shown in Fig. 123 The situation with the roles of Ve 
and reversed, i.e., v^. primarily connected to 1/3 and z/4, and primarily connected 
to Vx and z/2, gives similar results. Some examples of explicit 2 + 2 models are given 
in Refs. pimi]. 

The preceding arguments relied upon the assumption that large mixing solutions 
for solar or atmospheric neutrinos are satisfactory for oscillations to sterile neutrinos. 
In 2 + 2 scenarios, solar neutrinos oscillate to a linear combination of and i^^, and 
atmospheric neutrinos oscillate to the orthogonal combination |44j : 



— >■ — sin «!/,- + COS a Vg , (76) 

z/^ cos a z/t- + sin a Ug . (77) 

The oscillation probabilities in SuperK and SNO are |314j 

Rgj^ = pP + rpsin^ail- P), (78) 

RsNo = (3P, (79) 

RsNO = pP + Psm^a{l-P), (80) 



where R is the ratio of observed to expected rates in a given experiment, f3 is the 
ratio of the actual *B neutrino flux to the SSM prediction, P is the average oscillation 
probability for ^B neutrinos, and r = o'uf^u^/o'u^ — 1/6.48 is the ratio of z/^, i^r to 
z/e elastic scattering cross sections on electrons in the energy range of the SuperK 
experiment.'' The ratio of the SNO NC to CC data is then 

''SNO can also measure i/e scattering, which is equivalent to Rsk- 



72 



If the ®B neutrino flux is assumed to be known from the SSM [i.e., (3 = 1), then 
P is determined from SNO CC data and the sterile fraction can be deduced from 
the SNO NC/CC measurement. However, if P is allowed to be free, sin^ a cannot be 
determined since Rsno measures only the combination jSP and not P itself |314j . 

Early analyses of the solar neutrino data showed that solar solutions with pure 
^'e — > (a = 0) did not provide as good a fit to the solar neutrino data (see the second 
paper of Ref. [30] )• The main difficulty with pure sterile neutrino solutions is that 
they give similar predictions for the Chlorine and SuperK experiments, whereas active 
neutrino solutions give a larger value for SuperK due to neutral current i^^e and/or 

interactions in the detector, which is in better agreement with the experimental 
data. Oscillations to sterile neutrinos have different matter effects since the value of 
that gives resonant oscillations is = 2\/2GfEv{N(, — \Nn)/ cos 26*3 instead 
of the value 2\f2GFE^Ne/ cos29s for oscillations to active neutrinos (see Eg. I37|l. 

Imposing the SSM uncertainties on /3, the pure sterile solution is excluded at the 
7.6a level, although a large sterile fraction could still be allowed even after the SNO 
CC and NC data is included |107t I314| 1315] . The SuperK data (or, equivalently, the 
SNO ue scattering data) do not give an independent constraint; in fact, there is a 
sum rule [SHI Iffl 

RsNO = {RsK-il-r)R''sNo]/r, (82) 

that must be obeyed for any value of sin^ a. A future test of sin^ a that is independent 
of the SSM neutrino flux predictions would be a neutrino-nucleon NC measurement 
of intermediate energy solar neutrinos |314j . or the independent measurement of P, 
such as in the KamLAND reactor experiment |315j . 

The opposite extreme is to have pure sterile solutions to the atmospheric neutrino 
data {a = |). Here there are strong matter effects due to coherent forward scattering 
in the Earth that is present for z/^ but not Us (see Eq. |SZ|), with Neff = -\Nn; for 
pure —>■ oscillations, matter effects are small. The SuperK atmospheric data 
rule out pure Vs oscillations at the 99% C. L. |132[ 1317] . and require sin^ a to 

be smaller than 0.2 at the 90% C. L. |127| I131j . Other analyses also found that a 
substantial sterile component is allowed by the atmospheric neutrino data |318j . 

Since pure solar and pure z/^ — * Vg atmospheric neutrino oscillations are 
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ruled out, but partial sterile solutions are allowed in each case, the only remaining 
option in the 2 + 2 scenario is to have an active-sterile mixture for both solar and 
atmospheric neutrinos (0 < a < |). The analysis of Ref. jTZj indicates sin^ a > 0.55 
from pre-SNO salt phase solar neutrino data and sin^a < 0.35 from atmospheric 
neutrino data (both at 99% C. L.), and that the overall goodness of fit of 2 + 2 
models is only 1.6 x 10^^, much worse than that of the 3 + 1 models. 

To summarize the situation for four neutrinos, the most recent analysis jlTj shows 
that 3 + 1 scenarios provide a better fit than 2 + 2 scenarios, but neither scenario gives 
a good description of the combined solar, atmospheric, accelerator, and reactor data. 
However, the inclusion of additional sterile neutrinos can enhance the size of the LSND 
amphtude allowed by CDHS and Bugey in the 3 + 1 scenario ^Hl, and a recent study 
found that a 3 + 2 model does significantly better than the 3 + 1 models in fitting 
the data |319j . Also, studies using a more complete set of four-neutrino parameters 
indicate there may still be room for the 2 + 2 models j320j . Therefore, although 
the positive appearance result in LSND remains puzzling, a consistent four-neutrino 
explanation may still be possible. 

11.3 Three-neutrino models with CPT violation 

It has been suggested that if CPT were not conserved, then oscillations of three active 
neutrinos could describe the solar, atmospheric and LSND data simultaneously |321j . 
In this proposal, the mass matrices (and hence mass-squared differences) for neutri- 
nos and antineutrinos are different, which violates CPT (while preserving Lorentz 
invar iance) ^. 

In the original versions of CPT-violating models, the neutrino sector had the 

usual three-neutrino mass spectrum that can account for the oscillation of solar and 

atmospheric neutrinos, while in the antineutrino sector the mass-squared differences 

account for the oscillation of antineutrinos in the atmospheric and LSND experiments 

(the weak indication for z/^ — > z/g oscillations in LSND must be ignored). KamLAND 

data, consistent with oscillations of Ve at the 5m1 scale, forced a modification of 

^Whether such a model can be constructed using nonlocality of the interactions is stiU a matter 
of debate 
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the antineutrino spectrum, so that it describes the oscillation of antineutrinos in 
LSND and KamLAND (but not in the atmosphere) |323j . Since the atmospheric 
data does not distinguish between neutrinos and antineutrinos, this latter scenario 
was not in obvious contradiction to the data. However, one analysis of atmospheric 
data indicates that CPT- violating scenarios are not in good agreement with the 
atmospheric data |324j . Furthermore, global analyses of all data including KamLAND 
excludes these CPT- violating scenarios at the 3a level |324l I325j . 

We await the MiniBooNE results on z/^ oscillations to confirm or reject 

the LSND effect. A positive signal will rule out current CPT- violating models, and 
force a reconsideration of the disfavored four-neutrino models, possibly by extending 
them to include more than one sterile state. A negative result will rule out the 
standard four-neutrino scenario (which predicts the same oscillations for neutrinos 
and antineutrinos in LSND), but does not completely extinguish the LSND flame 
since it does not test the primary LSND —>■ z/g channel. In either case, it will be 
difficult to exclude the very speculative possibility of the LSND anomaly arising from 
CPT- violation in four-neutrino models |326j . 

12 Summary and outlook 

Great advances have been made over the past five years in our understanding of the 
neutrino sector. The field is now poised for further breakthroughs. In this article we 
have strived to summarize the present state of the field and the ongoing experimental, 
phenomenological, and theoretical efforts towards a bottom-up reconstruction of the 
fundamental properties and theory of neutrinos. In this concluding section we very 
briefly recapitulate the highlights of these recent accomplishments and the major 
planned directions for future progress. 

The standard formalism of three-neutrino mixing and oscillation probabilities was 
recounted. The modiflcations from Ue scattering in matter were quantifled as appro- 
priate for solar neutrinos and for long-baseline accelerator neutrino experiments. The 
experimental resolution of the solar neutrino problem, particularly by the SNO ex- 
periment, as the large mixing angle oscillation solution, was discussed along with its 
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validation by the KamLAND reactor and other neutrino experiments. The evidence 
for oscillations, from the dependence on the z/^ — z/^ oscillation probability 

versus zenith angle in the Super-Kamiokande atmospheric neutrino experiment, was 
summarized. 

In the framework of three-neutrino mixing, the major physics goal is to determine 
the six neutrino oscillation parameters: two mass-squared differences [Sm^, Sm^), 
three mixing angles {6s, Oa, 6x), and a CP phase 6. From present data the approxi- 
mate values are 6ml ^ 7 x 10"^ eV^, \6ml\ ^ 2 x IQ-^eV^, Os ^ 33°, ^ 45°. The 
unknown parameter 6x (with present limit 6x < 13° for 6ml = 2.0 x 10"^ eV^ at the 
95% C. L.) is critical to further understanding; it can be probed in z/^ or z/g z/^ 
appearance channels. Both off-axis and broad-band beams are being considered for its 
measurement. More intense neutrino beams (superbeams or neutrino factories) and 
longer baselines (> 1000 km) are essential to determine the sign of 6ml from matter 
effects and have sensitivity to CP-violation. Expected sensitivities of proposed future 
oscillation experiments to O^, sgn(5m^), and 6 are compared; see Table IHl Reactor 
experiments at baselines < 10 km could also determine Ox if sin^ 29x ^ 0.05, but only 
at about the 90% C. L. 

The evidence for — ^ z/g oscillations at a higher 6m^ (^0.1 — 1 eV^) may imply 
the existence of a sterile neutrino if CPT is conserved. Global fits to data now reject 
this possibility at more than the 99% C. L. The ongoing MiniBooNE experiment is 
designed to confirm or exclude the LSND effect. Big Bang Nucleosynthesis, combined 
with cosmic microwave background measurements by WMAP, strongly disfavors more 
than three thermalized neutrinos, but a large z/g-asymmetry allows an escape from 
this constraint. 

The search for the absolute scale of neutrino mass in beta decay, neutrinoless 
double-beta decay, and large scale structure in the universe was discussed. The best 
present constraint comes from cosmology and bounds the sum of neutrino masses to 
be less than about 1 eV. 

Galactic supernova explosions produce high fluxes of neutrinos and antineutrinos 
that can determine the neutrino mass hierarchy (normal or inverted) if the mixing 
angle Ox is known from accelerator experiments. 
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Models of neutrino mass are largely of two generic types, Grand Unification with a 
seesaw mechanism and radiative mass generation. Unified models can accommodate 
all present knowledge of quark and lepton masses and mixing, with different models 
giving different predictions for the unknown neutrino parameters. Successful Grand 
Unified models predict Sm^ > 0, for which there would be one heavier neutrino mass 
and two lighter masses. Also in unified models, the neutrino masses are hierarchi- 
cal: 1713 ^ "^2 ^ iTT-i- Radiative mass generation involves lepton number violating 
interactions that have testable implications for collider physics experiments. 

Unified models offer a possible explanation of the baryon asymmetry of the uni- 
verse in terms of a neutrino-antineutrino asymmetry. In such models the CP- violating 
phase in heavy Majorana neutrino decays that generate leptogenesis may be related 
to the CP- violating phase of neutrino oscillations at low energy. 

The ways by which the unknown neutrino oscillation parameters can be deter- 
mined by an ambitious future experimental program with baselines at long distances 
from intense sources was detailed. Conventional neutrino beam experiments are now 
under construction. Later these may be upgraded with more intense superbeams. 
Eventually neutrino factories will provide the ultimate sensitivity and precision in 
determining neutrino mixing and mass-squared difference parameters. 
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